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ABSTRACT 


In  the  study  of  fiber-reinforced  resin  composites,  the 
analysis  of  the  progressive  failure  of  a laminate  with  a stress  con- 
centration subjected  to  plane  stress  poses  a very  interesting  but 
complex  problem.  This  thesis  approaches  this  problem  by  using  the 
finite  element  method  to  examine  the  progressive  failure  of  symmetri- 
cal laminates. 

A modified  maximum  strain  failure  theory  is  proposed  and  a 
finite  element  computer  program  developed  that  accounts  for  progres- 
sive failure.  A computer  analysis  of  several  unnotched  laminate 
tensile  specimens,  with  lamina  at  various  angles,  was  made  and  these 
results  are  compared  with  experimental  data.  ^ 

, Circular  hole  tensile  specimens  with  (0°/90°/90°/0° ) and 
(0°/j_  45°/90° )^'  lamina  were  also  investigated,  and  the  progressive 
failure  through  the  finite  element  grid  presented.  The  ultimate 
failure  loads  of  the  circular  hole  specimens  are  cov\pared  with 
experimental  data.  Material  properties  used  were  those  for  Thornel 
300/5208  Graphite-Epoxy. 

Although  the  results  obtained  cannot  be  considered  conclusive 
for  all  cases,  they  do  compare  favorably  with  experiniental  data 
for  the  unnotched  specimens.  The  ultimate  failure  loads  of  ’■‘■j  hole 
specimens  were  somewhat  higher  than  those  obtained  experimentally. 
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INTRODUCTION 

The  word  "composite"  in  coi'iposite  material  signifies  that  two 
or  more  materials  are  combined  on  a macroscopic  scale  to  form  a use- 
ful material.  The  advantage  of  composites  is  that  the  materials  can 
be  combined  in  ways  that  usually  exhibit  the  best  qualities  of  their 
constituents  and  often  some  qualities  that  neither  constituent 
possesses.  Some  properties  that  may  be  improved  by  use  of  a compos- 
ite material  are  strength,  stiffness,  weight,  corrosion  resistance, 
fatigue  life,  and  thermal  properties.’ 

Composite  materials  have  been  used  for  centuries.  When  the 
first  composite  was  used  is  unknown,  but  recorded  history  contains 
references  to  various  forms  of  composite  materials.  For  example,  the 
Egyptians  used  laminated  wood  as  early  as  2780  B.C.,  and  the  Israel- 
ites added  chopped  straw  to  the  manufacture  of  bricks  in  800  B.C.^ 

A short  time  thereafter,  the  Mongol  bow  was  developed  from  a compos- 
ite of  animal  tendons,  wood  and  silk  bonded  together  with  an 
adhesive.  Still  later,  laminated  structures  appeared  in  the  Damascus 
gun  barrels  and  Japanese  ceremonial  swords.^ 

More  recently,  fiber-reinforced  resin  composites  that  have  a 
high  strength-to-weight  and  stiffness-to-weiqh‘  ratio  have  become 
important  in  weight  sensitive  applications  sui.n  as  aircraft  and  space 
vehicles.  Some  examples  of  these  modern  applications  of  fiber-resin 
composites  are:  an  AT-6C  aircraft  with  reinforced  plastic  fuselage 


2 


built  in  1943,'*  helicopter  rotor  blades  specifically  designed  to 
reduce  vibration  and  withstand  torsional  loads, rocket  motor  cases,' 
space  vehicle  structural  components,^  and  presently  an  increasing 
number  of  uses  such  as  fuselage  and  stabilizer  components  on  the  F-111, 
F-14,  F-15,  and  F-16  aircraft.'’®  Through  the  use  of  structural  com- 
pounds made  of  composite  materials,  strength-to-weight  ratios  have 
been  increased  100  percent  over  that  of  comparable  metal  structures.® 
The  impact  of  composites  on  jet  engine  performance  may  be  even  more 
dramatic,  where  an  800  percent  increase  in  the  thrust-to-weight  index 
appears  possible. ' 

The  superior  strength-to-weight  ratio  of  these  fibrous  composites 
is  related  to  the  failure  mechanisms  of  homogeneous  materials  where, 
generally,  the  actual  strength  is  considerably  lower  than  the  theo- 
retical atomic  strength.  The  reason  for  this  strength  difference  is 
the  formulation  and  movement  of  dislocations  in  the  homogeneous 
material.  By  forming  a material  into  thin  whiskers  or  fibers  with  a 
small  cross  section,  conformity  in  the  microstructure  is  enhanced,  the 
probability  of  internal  flaws  is  reduced,  and  the  formation  and 
movement  of  dislocations  restricted,  making  it  possible  for  the  fiber 
or  whisker  to  approach  its  theoretical  strength.®  Composite  sheets  or 
"lamina"  with  high  longitudinal  strength  are  formed  by  imbedding  many 
of  these  high  strength  fibers  longitudinally  in  a suitable  matrix 
material.  A composite  "laminate"  with  the  desired  strength  and  stiff- 
ness properties  may  then  be  formed  by  combining  layers  of  lamina 
together  at  various  orientations. 

Investigation  of  lamina  strennth  has  generally  been  approached 
from  both  micromechanical  and  macromechanical  levels. ’’ ® *' ® The 
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microinechanical  approach,  which  treats  a composite  material  as  a 
heterogeneous  continuum,  has  been  used  for  simple  lamina  model s ! ^ ^ ^ 
Although  theoretical ly  justifiable,  this  approach  has  a major  limita- 
tion in  that  the  analysis  required  is  extremely  complex,  and  there- 
fore, limited  to  very  simple  geometries.” 

In  general,  macromechanical  prediction  of  lamina  failure  has  been 
approached  from  one  of  the  following  three  theories:  the  maximum 

strain  theory,  maximum  stress  theory,  or  maximum  work  theory.'^”'*  Of 
these  three  theories,  the  maximum  work  approach  has  been  proven  to 
be  the  most  accurate  when  compared  with  experimental  data.'”^’'^ 
However,  this  theory  does  not  easily  lend  itself  to  the  analysis  of 
progressive  failure  because  the  damage  to  the  composite  cannot  be 
described  and  put  into  post-failure  relations.  The  maximum  strain  and 
maximum  stress  theories  are  well  suited  for  a progressive  failure  type 
analysis,  but  the  accuracy  of  these  theories  deteriorates  when  the 
fibers  are  at  an  angle  between  15  and  60  degrees  to  an  applied  uni- 
axial load.*  The  reason  for  this  loss  of  accuracy  at  intermediate 
angles  is  probably  due  to  not  considering  the  interactive  effect  on 
failure  of  combined  shear  and  tension. 

Failure  of  unnotched  laminates  has  been  approached  by  combining 
plies  through  lamination  theory  and  applying  a lamina  failure  theory 
to  each  ply.  In  doing  so,  the  disadvantages  of  ply  failure  theories 
are  carried  through  to  the  laminate.  In  addition,  failure  of  one 
lamina  as  a failure  criterion  for  the  laminate  is  usually  too  conser- 
vative. Maximum  work  or  distortional  energy  applications  to  laminate 
failure,  such  as  that  by  Tasi-Wu,^  have  given  good  results  for 
individual  laminates,  but  new  properties  must  be  obtained  for  each 
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1 

new  laminate.  Sendee kyj  usinq  the  method  of  Sandhu,^^  success- 
fully used  lamina  stress-strain  data  to  predict  the  nonlinear  response 
of  angle  and  multi-ply  laminates  in  uniaxial  tension.  The  failure 
‘ predictions  were  less  successful,  however.  The  success  of  this 

approach  to  progressive  failure  of  notched  laminates  is  still  to  be 
determined. 

• Because  of  the  complexity  of  the  micromechanical  approach, 

[ macromechanical  principles  are  usually  employed  to  determine  laminate 

i • 

behavior  in  the  presence  of  notches.  The  anisotropy  of  the  laminate 

f makes  the  failure  properties  due  to  stress  concentrations  of  particu- 

[ lar  interest  in  that  the  stress  concentration  factor  can  be  consider- 

! . I 

ably  higher  for  composites  than  for  isotropic  materials.  In  addition,  | 

a hole  and  crack  size  effect  has  been  observed  on  laminate  strength.  | 

Macromechanical  studies  of  notched  laminate  failure  have  used  models  ■ 

such  as  the  'inherent  flaw  model'  for  holes  by  Waddoups’^  and  an  | 

I 'average'  or  'point'  stress  approach  for  holes  and  cracks  by  Nuismer  i 

and  Whitney.^”  However,  both  studies  neglect  the  load-path  dependent  j 

i 

damage  or  progressive  failure  of  the  laminate.  In  doing  so,  none  can  ! 

I 

be  expected  to  have  general  applicability,  especially  when  biaxial  j 

I j 

I loading  is  considered.  j 

I The  purpose  of  this  thesis  is  to  study  progressive  failure  of  I 

I i 

notched  laminates  subjected  to  in  plane  loads.  i 

This  thesis:  j 

j 

(1)  presents  a modified  maximum  strain  theory  for  individual  ■ 

plies  and  develops  the  post  failure  constitutive  relations  j 

for  the  ply;  j 


(2)  develops  an  incremental  finite  element  program  which  uses 


i 


I 

( 

5 

laminated  plate  theory  to  account  for  stiffness  changes 
in  the  laminate  due  to  failure  in  the  plies; 

(3)  investigates  the  stress-strain  behavior  to  failure  of 

{ several  unnotched  laminates  (under  uniaxial  tension  loads) 

j and  compares  the  results  with  experimental  data; 

(4)  traces  the  progressive  failure  of  two  laminates  contain- 
ing circular  holes  and  compares  predicted  failure  loads 
with  experimental  data. 
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COMPOSITE  FAILURE 


Anisotropic  Elasticity 

With  the  advent  and  increased  usage  of  fibers  such  as  graphite 
in  composite  materials,  the  assumption  that  the  material  is  isotropic 
is  no  longer  valid.  Graphite  fibers  are  highly  anisotropic,  with  the 
longitudinal  stiffness  being  an  order  of  magnitude  greater  than  the 
transverse  stiffness.  Then,  in  order  to  analyze  such  fiber-rein- 
forced composites,  anisotropic  elasticity  must  be  employed. 

For  a three  dimensional  stress  state,  the  generalized  Hooke's 
Law  for  an  anisotropic  material  is  given  by: 


hi 

hi  ^ 

^2 

0-3 

• Wijl  < 

0-3 

V > 

^23 

^23 

^31 

"^31 

VI 2, 

V12, 

where  the  stiffness  matrix  is  symmetrical  with  21  independent 
constants. ' 

If  there  are  two  orthogonal  planes  of  material  property  symmetry, 
such  as  parallel  and  perpendicular  to  the  fibers  in  a unidirectional 
fiber  composite,  symmetry  will  also  exist  relative  to  a third  mutually 
orthogonal  plane,  and  the  material  is  said  to  be  orthotropic.  The  ^ 

stress-strain  relations  for  an  orthotropic  material  in  a coordinate 
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system  aligned  with  the  principal  material  directions,  or  parallel 
and  perpendicular  to  the  fiber  direction,  becomes:' 


''Oi  > 

Qll 

^12 

^13 

0 

0 

d~ 

^'0 

°2 

Qi2 

Q22 

^23 

0 

0 

0 

^2 

1^3 

^13 

^23 

^33 
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0 

0 

< 

^3 

1 ^23 

= 

0 

0 

0 

Q44 

0 

0 

^23 

^31 

0 

0 

0 

0 

Q55 

0 

^31 

0 

0 

0 

0 

0 

^66 

where  the  stiffness  matrix  has  been  reduced  to  nine  independent 
constants . 

Lamina  Constitutive  Relationships 

One  pertinent  assumption  in  establishing  the  constitutive  or 
stress-strain  relationships  for  the  lamina  of  a laminated  composite 
is  that  a lamina,  when  in  a composite,  is  in  a state  of  plane  stress. 
For  a state  of  plane  stress,  and  with  the  lamina  in  the  1-2  plane  as 
shown  in  Figure  1,  the  following  stresses  are  assumed  zero: 

03  = 0,  T23  = 0,  I.J3  = 0 (3) 
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Figure  1.  Unidirectional  Lamina 
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By  sutstitutinq  into  Equation  {2),  the  stress-strain  relation  becomes: 


where  the  components  of  the  stiffness  matrix  for  the  orthotropic  lamina 
qiven  in  terms  of  enqineerinq  constants  are: 


Q^l  - ^21^ 

^22  ~ ^22^^^”^!^  ^21^ 

Q-|2  ~ ^21^11^^^~^12  ^21^  ” ^12^22^^^”'^12  21^ 


(5) 


There  are  now  four  independent  constants:  E^-j,  E22’  ^’^2  ^12’ 

which  are  the  elastic  moduli  in  the  1 and  2 directions,  the  shear 
modulus  and  the  major  Poisson's  ratio,  respectively.  The  major  and 
minor  Poisson's  ratios  are  related  by: 

V2^Eii  = V.12E22  (6) 

where  the  major  Poisson's  ratio,  v.j2’  the  ratio  of  strain  in  the 
2 direction  to  strain  in  the  1 direction  for  a load  in  the  one  direc- 
tion. To  tailor  a material  with  the  proper  stiffness  and  strenqth 
in  various  directions,  unidirectional  laminae  are  usually  put  toqether 
with  fibers  runninq  in  several  different  directions  so  that  the  lamina 
principal  axes  are  not  coincident  with  the  reference  axes  of  the  lami- 
nate. When  this  occurs,  the  constitutive  relations  for  each  lamina 


I 

i 


must  be  transformed  to  the  laminate  reference  axes  in  order  to 


determine  the  laminate  constitutive  relationship.  The  transformation 
relationship  for  stress  between  an  arbitrary  x-y  axes  and  the  primary 
1-2  axes,  as  shown  in  Figure  2,  is 


while  that  for  strain  is 


/ \ 

/ \ 

^1 

e 

x 

^2 

"y 

l'^12 

1 

lY) 

where  the  transformation  matrix  is: 


[T,p  = 


-mn 


n 

2 

m 

mn 


2mn 

-2mn 

2 2 
m -n 


m = cos  0 
n = sin  0 


(7) 


(8) 


(9) 


axes  from  arbitrary  xy  axes 
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In  the  same  way,  the  primary  stress  and  strain  relations  are 
referenced  to  the  xy  axes  by: 


[T„J 


-1 


(10) 


(11) 


where  inverse  is  obtained  by  substituting  a negative  angle  G for 


the  positive  angle  0 in  the  T matrix.  Thus,  T inverse  becomes: 


[T,j] 


-1 


n 

mn 


- 2mn 
2mn 


2 2 

-mn  m -n 


(12) 


m = cos  0 
n = cos  0 


Knowing  the  orthotropic  lamina  material  properties  and  refer- 
encing them  to  the  xy  axes,  G is  measured  in  the  negative  direction. 
Then  T becomes  T'^  and  the  xy  stress-strain  relationship  obtained  from 
Equations  (4),  (8),  and  (10)  is 


-1 


= [T..]  [Q,^]  [T..]  c 


(13) 


where  Q-  is  the  orthotropic  lamina  stiffness  from  Equation  (4). 


11 


Denoting  the  lamina  stiffness  with  respect  to  tlie  xy  axes  as 


= [Q.J  e 


[Q. .]  = [T.  .]  [Q.  .]  [T.  .] 

‘■1.1  1J  1J  1.) 


Upon  multiplication,  the  terms  of  the  matrix  become: 

Qll  = + 2(Q^2  ^ 

Qi2  = (Q]]  Q22  ■ 

Qi3  = (Q^]  ^ ^12  ' ^^12  ' ^ZZ  ^ 

Q22  ^ ^ 2(Q^2  ^ ^22"'"^ 

^23  ""  ^"^11  '*'  ^12  ■ ^^12  ■ ^22 

^33  ^ ■*  ^22  ■ ^^^66^  ^66^"'^ 


A more  convenient  form  of  the  transformed  lamina  stiffness,  in  terms 
of  invariants  as  given  by  Tsai  and  Pagano^  and  used  later  in  the 
computer  analysis,  is: 

Qll  = U.|  + cos(2o)  + IJ^  cos(4o) 

^12  ""  ^4  ■ ^3 

Q'  .3  = i U,  sin(2n)  + U3  sin(4o) 

J3  2 2 3 

Q22  ' ■ *^2  ^3  cos(4o) 

Q23  ~ 2 ^Z  ■ *^3  sin(40) 

Q33  = Ug  - cos(4o) 
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where 


1 

' 8 

00,1 

+ 3Q22 

^ 2Qi2 

+ 4Q^ 

U2 

_ 1 
2 

(Q,, 

- Q22) 

^3 

1 

' 8 

(On 

+ Q22  - 

2Qi2  - 

^4 

1 

8 

(0,1 

+ Q22  + 

6Q12  - 

^5 

1 

' 8 

(0,1 

+ Q22 

2Qi2  + 

^^66 

Laminated  Plate  Theory 

Strain  Displacement  Relationshi p s 

A laminate  is  composed  of  several  orthotropic  layers.  As  such, 

the  description  of  the  behavior  of  a single  lamina,  as  previously 

discussed,  forms  the  basis  or  building  block  with  which  the  beh-’vior 

of  a laminate  may  be  described.  Equation  (14)  gives  the  constitutive 

relationship  for  a lamina  with  respect  to  an  arbitrary  xy  coordinate 

system.  Considering  the  arbitrary  xy  axes  to  be  oriented  with  the 

laminate  axes,  Equation  (14)  can  be  thought  of  as  a stress-strain 

t h 

relationship  for  the  k^  layer  of  a multi-layered  laminate  and  may  be 
written  as 


{o}|^  = [Q]|^{el|^  (19) 

Knowing  the  variation  of  stress  and  strain  through  the  laminate 
thickness  is  essential  to  the  definition  of  the  extensional  and  bend- 
ing stiffness  of  a laminate.  The  laminate  is  assumed  to  consist  of 
layers  of  perfectly  bonded  laminae,  such  that  the  displacements  arc 
continuous  across  lamina  boundaries  and  one  lamina  cannot  slip 
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relative  to  another.  With  this  assumption,  and  if  the  laminate  is 
thin,  it  may  be  assumed  that  a line  originally  straight  and  perpen- 
dicular to  the  middle  surface*  of  the  laminate  will  remain  straight 
and  perpendicular  to  the  inlodle  when  the  laminate  is  extended  and 
bent. 

Considering  a section  of  laminate  in  the  xy  plane  deformed  due 
to  some  loading,  as  shown  in  Figure  3,  the  geometrical  midplane  under- 
goes some  displacement,  u^,  in  the  x-direction.  With  the  above 
assumption,  the  line  ABD  remains  straight  and  normal  to  the  deformed 


undeformed 
cross  section 


Figure  3.  Geometry  of  deformation  in  the  xz  plane 


midplane  and  the  displacement  in  the  x-direction  of  any  point,  C,  on 
the  normal  ABD  is  given  by  the  linear  relationship^ 


Uc  = Uo  - 


(20) 
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where  is  the  z coordinate  of  the  ooint  C and  u is  the  slope  of  ABD 
with  respect  to  tlie  original  vertical  line.  Also,  under  deformation, 
line  ABD  remains  perpendicular  to  the  middle  surface  so  t' at  the  slope 
of  the  laminate  surface  in  the  x-di recti  on  is 


where  w is  the  displacement  in  the  z-direction.  Substituting  Equation 
(cl)  into  Equation  (20),  the  displacement,  u,  at  any  point,  z,  through 
the  laminate  thickness  is 


;iw 


By  similar  reasoning,  the  displacement,  v,  in  tiie  y-direction  is 

V = V - z (23) 

0 3y 

The  assumptions  thus  far  are  equivalent  to  ignoring  the  shearinn 
strains  in  planes  perpendicular  to  the  middle  surface,  that  is, 

= 0.  Also,  the  line  ABD  is  assumed  to  have  constant  length  so 
that  = 0.'  These  assumptions,  known  as  the  Ki rchhof f-Love  hypothe- 
sis, reduce  the  strain-displacement  relationships  to' 


hu  ^ Hv 
3y  Hx 


Differentiating  Equations  (22)  and  (23)  and  substituting  into  Equa- 
tion (24),  the  strains  become 


and  the  middle  surface  curvatures  are' 


By  substitution  of  the  strain  variation  throunh  the  thickness. 
Equation  (26)  into  Equation  (19),  the  stresses  in  the  k^*'  layer  can 
be  expressed  in  terms  of  the  laminate  middle  surface  strains  and 


m 


curvatures  as 
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^1.3 

/ 

X 

f'r  T 
X 

\ 

^22 

^23 

0 

^y 

> + 

y 

^'xy> 

1 

k 

^^13 

'^23 

^33 

k' 

1 

k 

0 1 

K 

xy/ 

/ 

]6 


(29) 


Laminate  Constituti Equa t i _qn_s 

Since  the  • can  be  different  for  each  layer  of  the  laminate, 
the  stress  variation  through  the  laminate  is  not  necessarily  linear 
even  though  the  strain  variation  is  linear.  To  investigate  these  non- 
linear stresses,  the  resultant  laminate  forces  and  moments,  denoted  by 
N and  M respectively,  are  obtained  by  integration  of  the  stresses  in 
each  layer  of  lamina  through  the  laminate  thickness.'  For  example. 


(30) 


where  t is  the  total  laminate  thickness. 

The  total  force  and  moment  resultants  for  a n-layered  laminate 
may  then  be  defined  as 


or,  using  Equation  (29)  and  summing  over  the  laminate  thickness,  for 
a laminate  with  n layers,' 


(33) 


(34) 


[c°]  and  [k]  are  not  functions  of  z,  but  are  midplane  values,  so  they 
can  be  removed  from  under  the  summation  signs  and  Equations  (33)  and 
(34)  can  be  written  as’ 


k=l 


I 

j 

I 

t For  laminates  that  are  symmetric  in  both  geometry  and  material 

i 

properties  about  the  middle  surface,  Equations  (35)  and  (36)  simplify 
considerably.  In  particular,  because  of  the  symmetry  of  [Q'^j]|/  ^nd 
the  lamina  thickness  t. , all  the  [B..]  are  equal  to  zero  and  the 

^ K 1 J 

force  and  moment  resultants  for  a symmetric  laminate  are 


For  the  remainder  of  this  investigation,  the  laminate  will  be 
considered  to  be  symmetrical  and  in  a state  of  pure  tension  or  com- 
pression, that  is,  bending  moments  will  be  zero  and  only  Equation  (40) 
will  apply. 
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PI y f.j  i 1 u rc-  C r it j'rjji 
Modified  Maximum  Stra i n TJioory 

Most  experimental  determinations  of  the  stremjth  of  a material 
are  based  on  uniaxial  stress  states.  However,  the  practical  problem 
usually  involves  at  least  a biaxial  state  of  stress.  For  an  ortho- 
tropic lamina,  strength  criteria  parallel  and  transverse  to  the  fiber 
direction  due  to  tension,  compression  and  shear  strength  may  all  be 
experimentally  determined.  To  relate  this  uniaxial  strength  informa- 
tion to  an  analysis  of  ply  damage  and  progressive  failure,  the  follow- 
ing modified  maximum  strain  theory  is  proposed.  The  lamina  is  said  to 
have  failed  in  the  fiber  direction  if 


^1  " ^t  ‘1  < ^ c 

and  transverse  to  the  fibers  if 


(42) 


'2  > Y^.t 


or 


to  < Y ^ 
2 cc 


(43) 


where  X X ,,  and  Y .,  Y ^ indicate  the  maximum  allowable  tensile 
ct  ec  > t cc 

and  compressive  strains  in  the  1 and  2 directions.  In  the  same  way, 
the  lamina  is  said  to  have  failed  in  shear  if 


'12 


(44) 


where  S is  the  maximum  allowable  shear  strain, 
c 

This  failure  theory  makes  it  possilile  to  obtain  post-failure 
constitutive  equations.  However,  stress  or  strain  interactions,  su'  h 
as  the  combined  effect  of  transverse  strain  and  shear  on  failure,  have 
been  ignored. 
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Post  Failure  Constituti v E cjuat j o n s 

In  forming  the  modified  strain  theory,  the  following  assumptions 
were  made: 

(1)  If  a lamina  fails  in  the  fiber  direction,  the  matrix  will 
still  carry  a load  transverse  to  the  fibers,  but  will  not 
carry  a shear  load. 

(2)  If  a lamina  fails  transverse  to  the  fiber  direction,  it 
will  not  carry  a shear  load,  but  the  fibers  will  carry  a 
normal  load  parallel  to  the  fibers. 

(3)  If  a lamina  fails  in  shear,  the  matrix  will  not  carry  a 
load  transverse  to  the  fibers,  but  the  fibers  will  carry 
a normal  load. 

(4)  If  a lamina  fails  in  the  fiber  direction  and  in  shear  or 
fails  both  parallel  and  transverse  to  the  fibers,  the  lamina 
is  considered  to  have  totally  failed  and  will  not  support 

a load. 

Each  of  the  above  assumptions  indicates  a partial  or  total  fail- 
ure of  the  lamina.  Examination  of  Equation  (4)  shows  that  for  a 
partial  failure  of  a lamina  at  a given  strain,  the  stress  is  changed 
by  a change  or  softening  of  the  stiffness  matrix.  In  the  computer 
solution  of  the  biaxial  stress  problem,  the  changes  in  stiffness  and 
load  are  used.  Therefore,  the  post-failure  lamina  constitutive  equa- 
tions will  be  expressed  in  terms  of  change  in  stress  and  stiffness  due 
to  partial  or  total  lamina  failure. 

Using  Equation  (4),  the  change  in  stress  due  to  a change  in 


stiffness  is  given  by 
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'ao^  ' 

Ao2 

‘2  > 

^12 
\ 7 

(45) 


where  is  the  change  in  stiffness  due  to  failure,  or 


AQ—  = - post  failure  stiffness 


(46) 


The  AQ^.j  terms  for  the  various  types  of  lamina  failure  are  found 
as  follows: 

(1)  Lamina  failure  in  the  fiber  direction  is  assumed  to  cause 
E^-|,  '^]2»  equal  zero  leaving  E22  as  the  only  factor 

contributing  to  the  new  stiffness.  Using  Equations  (5) 
and  (46),  the  AQ—  terms  are  found  to  be 


AQ 


11 


11  {~1  ” ^12'^^21^ 


AQ 


. ^22^'^12'^'2p 

22  (1  - 0^2^21  ^ 


AQ 


^’12  hi 


(47) 


12  rr^  ^’i2'^2i^ 


^^^66  ■ ^12 


(2)  In  the  same  way,  lamina  failure  transverse  to  the  fibers  is 
assumed  to  cause  622?  Q-|2'  ^-|2  equal  zero  leaving  E^^  as 
the  only  contributing  stiffness  factor.  For  this  type 
failure  the  AQ^j  terms  are 
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■ iTitM 


(1  - v-|2'^2p 


- '-'12^21^ 

^jaIu 

(i  - 


AQee  g^2 

(3)  Lainind  fdilure  in  shear  is  assumed  to  cause  £22*  "y^- 

to  equal  zero  leaving  E-j^  as  the  only  contributing  stiffness 
factor.  Thus,  Equation  (49)  also  gives  the  AQ^.j  terms  for 
shear  fai lure. 

(4)  Lamina  failure  in  the  fiber  direction  and  in  shear,  or  fail- 
ure both  parallel  and  transverse  to  the  fibers,  is  assumed 
to  cause  total  lamina  failure  and  therefore,  zero  remaining 
stiffness.  The  AQ^^  terms  obtained  by  Equations  (5) 
and  (46)  are 


11  (1  - v.|2V2i) 


^22  n - ^12^21 ^ 

. '’l_2_^ll_ 

^12  ■ (T  -“vY2'^‘2i'^' 


^^66  ' ^’12 


Denoting  AQ..  as  the  change  in  lamina  stiffness  with  respect  to 

* J 

the  arbitrary  xy  axes  and  using  Equations  (8),  (10),  and  (45),  the 
change  in  stress  in  the  xy  coordinate  system  due  to  a change  in  stiff- 
ness is  given  by 
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where 

[AQ,j]k  - (51) 

The  AQ^j  terms  are  calculated  by  use  of  Equations  (17)  and  (18)  where 
AQ--  terms  are  substituted  for  Q..  terms. 

^ vj  ' J 

Laminate  Failure  Criteria 

With  laminate  strength,  just  as  with  the  determination  of  lamin- 
ate stiffness,  the  basic  building  block  is  the  lamina  with  its  inherent 
characteristics.  Basic  to  determining  the  strength  of  a laminate  is 
a knowledge  of  the  stress  state  in  each  lamina.  However,  failure  of 
one  layer  does  not  necessarily  imply  failure  of  the  entire  laminate. 

The  laminate  may,  in  fact,  be  capable  of  higher  loads  despite  a 
significant  change  in  stiffness. 

The  strength  of  an  angle-ply  laminate,  symmetric  about  its  middle 
surface,  may  be  determined  by  examining  the  state  of  damage  in  each 
layer  for  a particular  load.  The  laminate  strains  are  calculated  from 
the  known  load  and  stiffness  prior  to  failure  of  a lamina.  If  one  or 
more  lamina  have  failed,  as  determined  from  the  failure  criterion,  a 
new  laminate  stiffness  is  calculated  and  the  laminate  strains  recalcu- 
lated to  determine  the  post-failure  strains.  Then  it  must  be  verified 
that  the  remaining  laminae,  at  their  increased  strain  levels,  do  not 


i 
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fail  at  this  applied  load.  Should  an  applied  load  cause  propressive 
failure,  where  all  layers  successfully  fail  at  the  same  load,  the 
laminate  is  said  to  have  suffered  gross  failure.* 

An  alternative  method,  described  in  the  next  section,  uses  the 
original  laminate  stiffness  to  determine  the  strains  at  each  load  or 
failure  cycle.  When  a failure  takes  place,  a change  in  stiffness  due 
to  the  failure  is  calculated.  Using  the  change  in  stiffness  and  the 
known  strains,  a pseudo  load  is  calculated  and  added  to  the  original 
load,  giving  the  required  increase  in  strain.  In  an  iterative  finite 
element  program  this  method  is  useful  in  that  the  stiffness  matrix  is 
only  inverted  once. 

Laminate  Post-Failure  Constitutive  Equations 

The  strength  of  a symmetric  angle-ply  laminate  subjected  to 
plane  stress  is  determined  by  first  finding  the  strains  for  a known 
load.  Inverting  the  stiffness  matrix.  Equation  (40)  can  be  written 
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From  the  previous  assumption  that  plane  sections  perpendicular 
to  the  midplane  axis  remain  plane,  and  for  a state  of  plane  stress. 
Equation  (52)  gives  the  state  of  strain  for  all  layers.  Then,  by 
Equation  (8),  the  strain  with  respect  to  the  1-2  axes  for  each  layer 
k,  is 


The  lamina  strains  are  compared  with  the  lamina  failure  criteria. 
Equations  (42),  (43),  and  (44),  to  determine  modes  of  failure. 


Should  failures  occur,  changes  in  stiffness  for  each  layer  are 

calculated  using  Equations  (47),  (48),  or  (49),  depending  on  the  type 

of  failure,  and  Equati""  (91).  The  total  change  in  laminate  stiffness 

is  found  by  summing  the  laminar  stiffness  changes.  For  an  n ply 

laminate,  the  total  change  in  stiffness,  '^E..,  is 

n 

[AE,j]  = 2 Wfjlk  '"k  - 

k=l 

where  AQ—  is  from  Equation  (51)  and  is  the  thickness  of 

lamina  n. 

Knowing  the  change  in  stiffness  and  the  laminate  strains,  a 
pseudo  force,  PN,  due  to  the  loss  of  stiffness  is  found  by 


Addim  this  pseudo  force  to  the  applied  load  of  Equation  (52)  gives 
the  increased  strain  due  to  lamina  failure.  Then,  the  new  strain 
is 
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(56) 


Equations  (53)  throuqh  (56)  are  repeated  until  equilibrium  is  obtained 
or  the  laminate  ex()eriences  qross  failure. 


Appl  ication  to  Coniposi  te^  Structures 

Thus  far  in  this  thesis,  lamina  and  laminate  have  been  analyzed 
in  a state  of  plane  stress,  but  the  qeometry  and  boundary  conditions 
have  not  been  considered.  It  has  been  assumed  that  the  state  of 
strain  is  constant  throughout  the  laminate  and  that  the  stress  is  con- 
stant throughout  each  layer.  In  actual  applications  the  stress  in  a 
laminate  and  in  a lamina  may  vary  considerably  due  to  qeometry  and 
loading  conditions.  Stress  concentrations  such  as  holes,  notches  and 
cracks  may  increase  the  local  stress  to  a much  greater  value  than  the 
stress  at  another  point  in  the  member.  Such  localized  stresses  may 
lead  to  localized  laminate  failure  and  ultimately  to  complete  laminate 
failure  at  reduced  loads. 

In  order  to  analyze  a varying  state  of  stress  at  points  across  a 


laminate,  the  finite  element  method  will  be  used  in  conjunction  with 
the  previous  ply  and  laminate  equations. 
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NUMERICAL  PROCEDURE 

FJ nj  te  El enient  Method  For  PJj_ne_  Stxesj  ^1^Ly§_L^_ 

In  a matrix  analysis  of  composite  materials  the  standard  approach 
is  to  divide  the  composite  laminate  into  a finite  number  of  elements 
connected  at  joints  or  nodal  points.  The  stiffness  or  flexibility 
properties  of  each  individual  element  are  then  established  by  an 
element  analysis,  and  the  element  stiffnesses  combined  to  form  the 
stiffness  matrix  for  the  complete  structure.  In  the  discussion  that 
follows,  a brief  description  of  the  displacement  method  for  a constant 
strain  triangle  element  will  be  presented  and  then  incorporated  with 
the  previous  ply  and  laminate  equations  in  an  iteration  method  to 
provide  a solution  to  the  nonlinear  composite  laminate  problem. 

Figure  4 depicts  a typical  triangular  element  with  nodes  i,  j, 
and  m,  numbered  in  counter-clockwise  order.  Each  node  may  have  dis- 
placements in  the  x and  y directions.  Then,  denoting  displacements 
in  the  x and  y directions  by  u and  v respectively,  the  six  components 
of  element  displacement  may  be  written  as  the  vector  {61  where 


'i  ■ “4  * '<5*1  * ‘Vi 


Vj  . t ..jKj  t 


V = a.  + iirX  + UcY 
111  4 5 m 6 m 


Evaluating  the  six  constants  a in  terms  of  the  nodal  displacements, 
gives^'  ^ 


u - 2^(a,+b,xtc,y)u.  * (Sj*b .xtc^yluj  * (a,„H,,x+c,„y)u„ 


where 


"i  = - Vi 


b ■ = y • - y 

1 j m 


S = %i  - " j 
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2jd=  det  i X-  y.  = 2 (area  of  triangle  ijm) 

J vJ 

^ ^m  -^m 


Neglecting  any  initial  strain,  the  total  strain  at  any  point 
within  the  element  can  be  defined  by  its  three  components  that  con- 
tribute to  the  internal  work.  From  Equations  (24) 
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Using  Equations  (57),  (61),  (62),  and  (65),  the  strain  within  the 
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Nodal  forces  can  also  be  expressed  in  terms  of  the  nodal  dis- 
placements. Denoting  U and  V as  the  nodal  forces  in  the  x and  y 
directions  and  assuming  zero  body  forces,  the  nodal  force  vector  iF) 


for  a triangular  element  can  be  written 
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(70) 


The  stresses  that  result  from  these  nodal  forces  can  be  found 
by  equating  the  work  done  by  the  forces  to  the  strain  energy  stored  in 
the  element.  The  work  done  by  the  nodal  forces  is^^ 


u = 1 (U,.,  + V,»,  * U.Uj  + VjVj  * * V„vJ 


or 


W = \ (Fl® 


The  strain  energy  is  given  b/ 


II  = — r f (a  c + OC  + i Y )dA 
2 J J ' X X y y xy'xy' 
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(71) 


(72) 


(73) 


where  .d  is  the  area  of  the  triangular  element  and  t is  the  element 
thickness.  Using  Equations  (66)  and  (69),  the  strain  energy  can  be 
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rewri tten 

U = [B]’’’[E]^[B]{6}®  (74) 

Equating  the  work  and  energy  equations  and  taking  the  transpose  of 
both  sides, 

{F}®  =^t  [B]^[E][B]{6}®  (75) 

Denoting  the  element  stiffness  matrix  [K]^, 

{F}®  = [Kj®{6}^  (76) 

and 

[K]®  =Zlt  [B]^[E][B]  (77) 

Equations  (76)  and  (77)  are  now  sufficient  for  computation  with  the 
actual  matrix  operations  being  accomplished  in  the  computer  program. 
Combining  the  element  stiffness  matrices  and  their  force  and  dis- 
placement vectors  gives  the  structural  system  of  equations 

[A]  {6}  = {F}  (78) 

where  [A]  is  the  structural  stiffness  matrix. 

Solution  Method  foi^  Nonlinear  Material  Properties 
Initial  Stress  Process 

The  expressions  derived  in  the  previous  sections  describe  fully 
the  stress-strain  relations  for  a laminated  composite  material  in  a 
state  of  plane  stress.  The  essential  nonlinearity  is  evident  from 
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Equations  (54),  (55),  and  (56)  with  the  composite  stiffness  matrix 
being  dependent  on  the  state  of  total  stress.  This  problem,  as 
described  in  the  following  section,  can  be  approached  using  peace- 
wise  linearization  to  obtain  a solution  iteratively.'”’^® 

The  "initial  stress"  process  approaches  the  solution  of  a non- 
linear problem  as  a series  of  approximations . ‘ ^ ® In  the  first 
step  after  a load  increment  a purely  elastic  problem  is  solved  deter- 
mining an  increment  of  strain  {Ac'}  and  of  stress  {Ao' } at  every  point 
of  the  structure.  The  nonlinearity  implies  that  for  the  increment  of 
strain  found,  the  increment  of  stress  will,  in  general,  not  be 
correct.  If  the  true  increment  of  stress  for  equilibrium  is  {A'\l,  then 
the  correct  solution  can  be  maintained  by  a set  of  pseudo  body  forces 
equilibrating  the  "initial  stress"  system  (Ao'l  - fAa}.'” 

At  the  second  stage  of  computation  the  system  of  pseudo  body 
forces  can  be  removed  by  allowing  the  structure  (with  unchanged  elastic 
properties)  to  deform  further.  An  additional  set  of  strain  and  stress 
increments  is  caused,  and  once  again  they  are  likely  to  exceed  those 
permitted  by  the  nonlinear  problem.  The  redistribution  of  pseudo  body 
forces  is  repeated  and  the  process  continued  until  it  converges  to  the 
nonlinear  equilibrium  conditions. 

Application  to  Composite  Materials 

In  laminated  composite  materials,  the  nonlinearity  comes  from 
failure  or  partial  failure  of  a ply  within  a laminate.  Ply  failure 
or  partial  failure  implies  that  a change  in  stiffness  has  taken  place 
and  that  the  load  used  and  displacements  found,  for  an  elastic  solu- 
tion, are  not  correct.  To  arrive  at  the  correct  solution,  pseudo  body 
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forces  are  calculated  using  the  change  in  stiffness  and  the  laminate 
strains.  These  pseudo  forces  are  allowed  to  further  deform  the  lami- 
nate using  the  original  elastic  properties.  New  strains  are  found  and 
the  process  repeated  until  equilibrium  is  obtained. 

Specific  steps  in  the  initial  stress  process  as  applied  to 
composite  materials  are; 

(1)  The  problem  is  set  up  by  using  Equation  (77),  for  each 
element,  to  construct  the  structural  stiffness  matrix. 

An  incremental  load  and  other  boundary  conditions  are 
entered  into  Equation  (73)  which  gives  the  following 
system  of  equations  to  be  solved: 

[A]  {cS}  = {FI  (79) 

(2)  The  stiffness  matrix  [A]  is  partially  inverted  and  the 
displacement  {6}  computed. 

(3)  Strains  within  each  element  are  found  by  Equation  (66) 


(4)  From  Equation  (8),  the  principal  strains  in  each  lamina 
are  obtained, 
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(5)  Lamina  failure  criteria.  Equations  (42),  (43),  and  (44) 
are  applied.  If  there  are  no  failures,  qo  to  step  10. 

(6)  The  change  in  element  stiffness  due  to  failure  is  computed 

using  Equation  (54), 
n 

[AE]  = [AQ]^,  (z^  - z^_^)  (82) 

k=l 


(7)  Pseudo  forces  at  each  node  point  are  found  using  [aE]  as 
the  element  stiffness  and  Equations  (66)  and  (75).  Denot- 
ing the  pseudo  forces  PF, 


(PFl  = [B]^[AE] 


(83) 


(8)  The  maximum  pseudo  force,  is  compared  with  an 

accuracy  constant  ACC.  If  PF^^^^  is  less  than  ACC,  equilib- 
rium has  been  reached,  go  to  step  10. 

(9)  Using  the  pseudo  forces  and  the  partially  inverted  stiff- 
ness matrix  of  step  2,  additional  displacements  are  found 
and  added  to  the  original  displacements.  Return  to  step  3. 

(10)  The  load  is  incremented  by  adding  the  displacements 

obtained  for  an  incremental  load  to  existing  displacements 
and  returning  to  step  3. 

Should  the  iteration  process  of  steps  3 through  9 be  repeated 
20  times  within  an  increment  without  reaching  equilibrium,  the  laminate 
is  considered  to  have  suffered  gross  failure  and  the  process  is 
stopped. 
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SOLUTION  OF  PROBLEMS 


Uniaxial  Tension  SjJGcinienj^ 

To  check  the  reliability  of  the  modified  strain  theory  and  the 
finite  element  program,  the  predicted  stress-strain  curves  and  failure 
loads  for  (0”)^,  (90")^,  (0'790V90V0'’ )j. , (90°/+  45^0°  and 
(90°/+  45°/90° laminates  were  obtained  for  uniaxial  tension  loads. 
Figure  5 shows  the  finite  element  grid  used.  The  12  element  grid, 
scaled  to  5 inches  (12.7  cm)  in  length  and  1 inch  (2.54  cm)  in  width, 
was  loaded  using  incremental  displacements  in  the  direction  shown. 

Zero  displacement  conditions  were  specified  for  nodes  opposite  the  load 
end  in  the  load  direction  and  along  one  side  transverse  to  the  load 
direction. 

To  establish  a stress-strain  relationship  for  comparison  to 
experimental  data,  one  element  was  chosen  and  its  state  of  stress  and 
strain  written  out  at  the  end  of  each  increment.  The  failure  status 
of  eacy  ply  within  each  element,  nodal  displacements,  iterations  and 
the  maximum  pseudo  force  for  each  iteration  were  also  written  out. 

Material  properties  used  were  those  for  Thornel  300/5208  graphite- 
epoxy,  listed  in  Appendix  A. 

yjicul  a_r  ]lql  e_  Spec  ntjens 

Two  laminates  containing  circular  holes  and  loaded  in  uniaxial 


tension  were  investigated  using  the  finite  element  arid  shown  in 
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Figure  6.  The  99  element  grid  was  given  an  incremental  displacement 
load  in  the  direction  shown,  with  zero  displacement  conditions  imposed 
on  the  end  opposite  the  load  in  the  load  direction  and  along  the  hole 
side  transverse  to  the  load  direction.  The  scale  for  the  grid  repre- 
sented a specimen  3 inches  (7.62  cm)  wide  with  a hole  1 inch  (2.54  cm) 
in  diameter. 

At  the  end  of  each  increment  the  status  of  eacy  ply  within  each 
element  was  printed  out.  Total  failure  loads  were  determined  using 
the  nodal  displacements  and  stiffnesses  at  failure  to  calculate  the 
nodal  forces  at  the  load  points. 

Material  properties  used  were  those  for  Thornel  300/5208  graphite- 
epoxy,  listed  in  Appendix  A. 
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RESULTS 


Uniaxial  Tension  Specimens 

Stress  vs.  strain  diagrams  for  the  uniaxial  tension  problems 
are  given  in  Figures  7 through  11.  Data  points  plotted  are  those 
obtained  by  Sendeckyj  for  Thornell  300/5208  graphite-epoxy  laminates’" 
The  solid  line  is  the  predicted  stress-strain  curve  as  obtained  from 
a chosen  element.  Changes  in  nodal  displacements  written  out  at  the 
end  of  each  increment  v/ere  egual  for  all  elements,  indicating  the 
strains,  and  thus  the  stresses,  were  equal  for  all  elements. 

c jji'  e ti  s 

(0V90V90°/0'’)^  Laminate 

Figures  12  through  16  give  the  damage  or  progressive  failure 
status  of  the  (0°/90"/90°/0°)^  laminate  atthe  end  of  each  load  incre- 
ment. The  first  number  code  indicates  the  status  of  the  O'’  plies 
and  the  second  number  the  status  of  the  90"  plies.  The  specific 
numbers  give  the  mode  of  failure  within  the  ply. 


The  total  failure  load  calculated  for  this  notched  laminate  was 
37,600  psi  (2.6  x 10^  Pa)  while  that  obtained  experimentally  by 
Nuismer  and  Whitney  was  28,200  psi  (1.9  x 10^  Pa).’“ 


strain  (10  ) 


Figure  10.  Stress  vs  strain,  (90"/+  45'’/0°)^  laminate  in  uniaxial 

tension 


stress 


0 = no  failure  2 = matrix  failure 

1 = fiber  failure  3 = matrix  and  fiber  failure 


Figure  12.  (0'’/90°/90'’/0'’ circular  hole  partial  failure  at 

.009  in.  (.023  cm)  displacement  load 


1 

51 

(0°/l  45790°)^  Laminate 

Figures  17  through  20  give  the  damage  or  progressive  failure 
status  of  the  (0°/+  45°/90° laminate  at  the  end  of  each  load  incre- 
ment. The  numbers  in  the  code  indicate  the  failure  status  of  the  plies 
as  follows;  the  first  number  for  the  0”  ply,  the  second  number  for 
the  + 45°  ply,  the  third  number  for  the  - 45°  ply  and  the  fourth 
number  for  the  90°  ply. 

The  total  failure  load  was  calculated  to  be  47,000  psi  (3.3  x 
10^  Pa)  while  an  experimental  value  of  45,700  psi  (3.2  x 10^  Pa) 
was  obtained  by  Nuismer  and  Whitney.’” 
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0 = no  failure  2 = matrix  failure 

1 = fiber  failure  3 = matrix  and  fiber  failure 


Figure  18.  (0”/+  45”/90” circular  hole  partial  failure  at 

.012  in.  (.030  cm)  displacement  load 
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0 = no  failure  2 = matrix  failure 

1 = filjer  failure  3 = matrix  and  fiber  failure 


Figure  19.  (0°/+  45”/90° circular  hole  partial  failure  at 

.016  in.  (.041  cm)  displacement  load 


0 = no  failure  2 = matrix  failure 

1 = fiber  failure  3 = matrix  and  fiber  failure 


Figure  20.  (0’’/+  45°/90'’ circular  hole,  complete  failure  at 

.020  in.  (.051  cm)  displacement  load 


DISCUSSION 


In  qeneral , the  theoretical  stress-strai n results.  Figures  7 
through  11,  for  the  unttotcheci  tensile  specimens  compare  favorably 
with  experimental  data  with  some  variation  in  the  ultimate  failure 
point  of  the  laminate.  Figures  7 and  8 show  the  theoretical  failure 
stresses  and  strains  for  the  0°  and  90°  laminates  to  be  somewhat  below 
the  experimental  failure  values.  These  differences  might  be  the  result 
of  choosing  low  values  for  the  maximum  allowable  strains.  However, 
Figures  9 and  11,  using  the  same  maximum  allowable  strains,  indicate 
theoretical  ultimate  strengths  for  the  (0°/90°/90°/0° and 
(90°/+  45°/90°)g  laminates  slightly  higlier  than  the  experimental 
values.  The  low  experimental  failure  values  of  Figure  9 were  sus- 
pected to  be  due  to  damage  to  the  surface  0°  plies  during  specimen 
handling  and  fabrication.  In  Figure  11,  where  the  +^45°  plies  are 
the  dominant  load  carrying  plies,  deviation  from  the  experimental 
values  may  be  influenced  by  the  assumption  that  the  shear  stress  vs 
strain  is  linear,  when  in  actuality,  it  is  highly  nonlinear. 

(90°/+  45°/90° theoretical  results.  Figure  10,  compare  well  with 
experimental  data.  The  horizontal  jumps  in  the  theoretical  curves  of 
Figures  9,  10,  and  11,  indicate  increased  strain  due  to  failure  of 
the  90°  plies. 

The  progressive  failure  of  the  circular  hole  specimens  produced 
interesting  results  in  that  the  failure  of  both  specimens  beqan  at 
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the  hole  edge  in  a direction  perpendicular  to  the  load  direction  but  ‘j 

did  not  progress  in  the  shortest  direction  to  the  specimen  outer  edge.  i 

Figures  li^  througli  16,  and  Figures  17  througli  20,  seem  to  indicate 
that  for  both  the  (0"/90"/90V0°)^  and  (0°/i  45790° laminates 
the  failure  path  was  approximately  45  above  the  shortest,  or  hori-  I 

zontal  path  across  the  grid.  However,  it  should  be  noted  that  although 
the  failure  status  of  the  elements  in  Figures  16  and  20  probably  give 
a good  indication  of  the  ultimate  failure  modes,  they  may  not  be  exact 
because  equilibrium  of  the  pseudo  nodal  forces  was  not  obtained  for  | 

the  failure  load. 

Figures  12  through  16  show  the  90°  plies  of  the  (0°/90°/90°/0° 
laminate  failing  first  in  the  matrix  as  expected,  with  the  0°  plies 
following  essentially  the  same  element  failure  pattern  at  higher  loads. 

The  ultimate  failure  load  of  37,600  psi  (2.6  x lO''^  Pa)  obtained  theo- 

g 

retically  was  much  higher  than  the  28,200  psi  (1.9  x 10  Pa)  experi- 
mental value.  This  error  of  over  30  percent  possibly  indicates  that 
the  modified  strain  theory  used  needs  refinement. 

Figures  17  through  20,  for  the  (0°/j^  45°/90° laminate  show 
the  90°  ply  failing  first  in  the  matrix,  with  this  failure  oroqressino 
throughout  most  of  the  structure  before  the  ultimate  failure  load  was 
reached.  Failure  of  the  + 45°  and  0°  plies  is  more  limited  throughout 
the  load  range  and  indicates  the  progressive  path  of  total  failure. 

Although  the  theoretical  failure  load  of  47,800  psi  (3.3  x 10^  Pa)  is 

O 

close  to  the  experimental  value  of  46,700  psi  (3.2  x 10'"  Pa),  this 
experimental  value  is  suspected  to  be  too  high.  The  error  is  sus- 
pected because,  in  the  study  by  Nuismer  and  Whitney,^” 


other  failure 


stresses  decreased  with  incroasincj  hole  size  as  expected,  whereas 
this  particular  value  increased.  Comparison  with  other  data  in  this 
study^°  indicates  that  the  actual  failure  stress  miqht  be  around 
40,000  psi  (2.8  X 10“^  Pa).  If  this  is  the  case,  the  error  is  consider- 
able and  perhaps  is  again  pointing  to  a necessary  refinement  of  the 
modified  strain  theory. 

In  all  problem  solutions,  it  v/as  assumed  that  the  lamina  had  the 
same  stress-strain  curve  in  compression  as  tension.  Also,  compressive- 
failure  strains  were  assumed  to  be  the  same  as  those  for  tension. 

These  assumptions,  made  partially  because  of  the  lack  of  reliable  com- 
pression data,  are  certainly  incorrect  and  would  have  to  be  modified 
for  problems  involving  substantial  compression.  However,  for  the  hole 

I 

I 

I problems  considered,  where  compression  occurs  only  in  a small  region 

at  the  top  of  the  hole,  this  is  not  expected  to  result  in  the  appre- 
ciable errors. 

One  other  factor  not  considered  in  the  investigation  is  the  free 
edge  effects. The  assumption  was  made  that  the  strain  through  the 
laminate  was  constant  at  any  given  point.  With  this  assumption,  the 
stresses  in  plies  at  different  orientations  will  generally  be  differ- 
ent. At  free  edges,  although  the  average  stress  along  the  free  edge 
is  set  equal  to  zero,  this  leads  to  mathematically  non-zero  surface 
tractions  along  the  free  edge  of  each  ply,  thus  violating  the  actual 
boundary  conditions.  If  the  actual  boundary  conditions  are  used,  it 
can  be  shown  to  result  in  significant  interlaminar  shear  and  normal 
stresses  that  have  been  shown  to  be  responsible  for  delamination  along 


free  boundaries.’^  However,  this  effect  has  been  shown  to  perturb 
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the  inplane  stresses  predicted  from  laminated  plate  theory  in  only 
a small  region  near  the  boundary.  Furthermore,  observation  of  notched 
specimens  subjected  to  monotonic  failure  loads  has  produced  no  evi- 
dence of  delamination  at  the  notch  before  failure  occurs.  Thus,  the 
free  edge  effect  is  expected  to  be  of  importance  in  the  progressive 
failure  of  notched  laminates  only  if  fatigue  loadings  are  considered. 
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COMPUTER  PROGRAM  DESCRIPTION 
Maj  n_  P ram 

The  Main  program  is  the  executive  routine  that  controls  the 
sequence  of  steps  by  cal  lino  subroutines  that  set  up  and  execute 
the  problem.  Subroutines  called  by  Main  are  given  in  the  Computer 
Program  section. 


Subroutines 


Setup 

Subroutine  Setup  is  called  by  the  Main  program.  It  reads  the 
data  deck,  adjusts  the  structure  size  by  scaling  factors,  checks  the 
bandwidth  and  adjusts  it  if  necessary.  Setup  then  calls  Stifgn,  and 
upon  return  of  control,  writes  out  the  data  deck  and  returns  control 
to  Main. 

Stifgn 

Subroutine  Stifgn  computes  the  laminate  stiffness  matrix  by 
using  the  orthotrooic  lamina  properties  and  Equations  (17),  (18),  and 
(37).  Control  is  returned  to  Setup. 

Constr 

Subroutine  Constr.  is  called  by  Main.  It  controls  the  construc- 
tion of  the  system  stiffness  matrix  by  calling  Subroutine  Elcons  for 
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each  element,  and  then  returns  control  to  Main. 

Elcons 

Subroutine  Elcons  is  called  by  Constr.  This  subroutine  con- 
structs the  stiffness  matrixes  for  the  individual  elements  by  usinq 
Equation  (77),  and  then  combines  the  element  stiffnesses  to  form  the 
structure  stiffness  matrix.  Control  is  returned  to  Constr. 

ExcUe 

Subroutine  Excite  is  called  by  the  Main  proqram  and  enters  tiie 
problem  boundary  conditions.  If  the  boundary  conditions  are  specified 
forces,  these  forces  are  entered  directly  into  the  force  matrix.  For 
displacement  boundary  conditions,  a pseudo  force  dependent  only  on 
the  specified  displacement  is  entered  into  the  force  matrix,  the 
corresponding  diagonal  elements  of  the  stiffness  matrix  are  not 
changed,  but  assumed  equal  to  one,  and  the  remaining  terms  in  the 
related  rows  and  columns  are  set  equal  to  zero.  Control  is  returned 
to  Main. 

Gausel 

Subroutine  Gausel  is  called  by  the  Main  program.  It  solves  the 
boundary  value  problem  by  Gaussian  elimination,  leaving  the  stiffness 
matrix  in  partially  inverted  form.  This  partially  inverted  matrix  is 
subsequently  used  in  Subroutine  Foredu  during  iteration.  Control  is 


returned  to  Main. 
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EJ/ail 

Subroutine  Elfail  is  called  by  the  Main  program.  This  sub- 
routine controls  the  system  failure  analysis,  iteration,  incremen- 
tation and  output.  Subroutines  called  by  Elfail  are:  Strain,  Layer, 

Pforce,  Foredu,  and  Output.  Upon  failure  to  reach  equilibrium 
within  20  iterations,  control  is  returned  to  Main. 

Strain 

Subroutine  Strain  is  called  by  Elfail.  It  calculates  element 
strains  using  Equatiofi  (66)  and  returns  control  to  Elfail. 

Layer 

Subroutine  Layer  is  called  by  Elfail.  This  subroutine  applies 
the  failure  criteria  to  each  lamina  within  each  element.  Should 
failure  occur,  it  is  noted  in  the  failure  tracing  array  ITT,  Subrou- 
tine Dlstif  is  called,  then  control  is  returned  to  Elfail.  If  no 
failures  occur,  control  is  returned  to  Elfail. 

Dlstif 

Subroutine  Dlstif  is  called  by  Layer.  It  calculates  the  change 
in  element  stiffness  due  to  failures  determined  in  Layer.  Co‘'cula- 
tions  are  made  using  Equation  (54)  and  control  returned  to  Layer. 

Pforce 

Subroutine  Pforce  is  called  by  Elfail.  Using  Equation  (83),  this 
subroutine  calculates  the  pseudo  nodal  forces  due  to  changes  in  stiff- 
ness. Control  is  returned  to  Elfail. 


Subroutine  Foreciu  is  called  by  El  fail, 
inverted  stiffness  matrix  to  reduce  the  force 
substitutes  to  solve  for  displacements  due  to 
tation.  These  displacements  are  added  to  the 
and  control  is  returned  to  Elfail. 


It  uses  the  partial ly 
matrix  and  then  bacE 
iteration  or  incremen- 
total  displacement  vector 


Output 

Subroutine  Output  is  called  by  Elfail.  At  the  end  of  each 
increment  or  upon  failure  to  reach  ec|ui  1 i bri urn.  Output  writes  out 
the  applied  displacement  or  load,  the  failure  status  of  each  lamina 
within  each  element,  the  nodal  displacements  and  the  pseudo  nodal 
forces.  Control  is  returned  to  Elfail. 
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APPENDIX  C 

COMPUTER  PROGRAM  VARIABLES 

A System  stiffness  matrix 

FO  Force  matrix,  nodal  values 

FIN  nitial  force  matrix  with  boundary  conditions  entered, 

used  in  incrementation 

DIS  Delta  displacement  vector,  nodal  values 

TDIS  Total  displacement  vector,  nodal  values 

Q Orthotropic  lamina  stiffness 

EE  Delta  stiffness  due  to  failure 

Ell  Ell  through  E33  make  up  the  laminate  stiffness  matrix 

E13 

E22 

E23 

E33 

ANG2R  Two  times  lamina  angle  in  radians 

THETA  Lamina  angle  in  radians 

ITT  Failure  tracing  matrix 

STRN  Strain  in  element  chosen  for  output 

STRS  Stress  in  element  chosen  for  output 

FOO  Force  matrix  use  in  iteration 


FOMAX  Maximum  nodal  force  during  iteration 
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EPX  Element  strain  in  the  x direction 

EPY  Element  strain  in  the  y direction 

GXY  Element  shear  strain 

ESTRN  Lamina  strain  matrix  referenced  to  xy  coordinate  system 

LMSTRN  Lamina  strain  matrix  referenced  to  the  principal  lamina  axes 

LMSTRS  Lamina  stress  matrix  referenced  to  the  principal  lamina  axes 


B 

B transpose 

used  in 

Subroutine  Pforce 

BB 

Calculation 

matrix 

used 

in 

Subroutine 

Pforce 

EP 

Calculation 

matri x 

used 

i n 

Subroutine 

Pforce 

APPENDIX  D 


IBD 

NRD 

NEL 

NXY 


XYM 

SEX 

SFY 

NXD 

NXF 

NXl,  2 
NSX 

NYD 

NYF 

NYl,  2 
NSY 

NZC 

NANG 


COMPUTER  PROGRAM  INPUT  DATA 

Bandwidth;  2(largest  difference  in  node  numbers  + 1) 

Matrix  order;  2(number  of  nodes) 

Number  of  elements 

1;  Coordinate  data  is  input  sequentially 
2;  Coordinate  data  is  input  as  one  pair  per  card 

Nodal  numbers  for  each  element  in  counter-clockwise  order 

Nodal  coordinates,  ordered  pairs,  x-coordinate  first 

X-scalinq  factor 

Y-scaling  factor 

Number  of  non-zero  applied  displacements,  x-face 

Number  of  non-zero  applied  forces,  x-face 

End  points  of  integration  path  to  get  a total  force,  x-face 

1;  Resulting  force  is  based  on  loads  applied  to  the  x-face 

2;  Resulting  force  is  based  on  displacements  applied  to 

the  x-face 

Number  of  non-zero  applied  displacements  y-face 

Number  of  non-zero  applied  forces,  y-face 

End  points  of  integration  path  to  get  a total  force,  y-face 

1;  Resulting  force  is  based  on  loads  applied  to  the  y-face 

2;  Resulting  force  is  based  on  displacements  applied  to 

the  y-face 

Total  number  of  zero  displacements 
Number  of  unique  ply  orientations 


NFAIL 


1 ; Maximum  strain  failure 
2;  Maximum  stress  failure 
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NDPX 

NFPX 

EX  (1)  = 

EX  (2)  = 

NX 

NDPY 

NFPY 

EY  (1)  = 

EY  (2)  = 

NY 
NZP 

El 

E2 

G 

V12 

ANGLE 

THK 

ALLOWd  ,1  ) = 
ALLOW(2,l )= 
ALL0W(3,1 )= 
ALL0W(4,1 )= 


Array  positions  of  coordinate  numbers  of  non-zero  applied 
displacements,  x-face  (2  X node  number  -1) 

Array  positions  of  coordinate  numbers  of  non-zero  applied 
forces,  x-face  (2  X node  number  -1) 

Maepiitudo  of  applied  displacement  increment,  x-face 

Magnitude  of  applied  force  increment,  x-face 

Nodal  numbers  adjacent  to  applied  force  nodes,  x-face 

Array  positions  of  coordinate  numbers  of  non-zero  applied 
displacements,  y-face  (2  Y node  number) 

Array  positions  of  coordinate  numbers  of  non-zero  applied 
forces,  y-face  (2  X node  number) 

Magnitude  of  applied  displacement  increment,  y-face 

Magnitude  of  applied  force  increment,  y-face 

Nodal  numbers  adjacent  to  applied  force  nodes,  y-face 

Array  identification  numbers  for  zero  displacement 
condi ti ons 

Orthotropic  material  modulus  in  fiber  direction 

Orthotropic  material  modulus  transverse  to  fibers 

Orthotropic  material  shear  modulus 

Orthotropic  material  major  Poisson's  ratio 

Orientation  angles  of  individual  plies,  positive  counter- 
clockwise from  the  x-axis,  in  degrees 

Thickness  of  all  plies  at  each  unique  orientation 

Limiting  ply  tensile  strain,  parallel  to  fibers 

Limiting  ply  tensile  strain,  transverse  to  fibers 

Limiting  ply  shear  strain 

Limiting  ply  tensile  stress,  parallel  to  fibers 


1 
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ALL0W(5,1)=  Liinitinq  ply  tensile  stress,  transverse  to  fibers 
ALL0W(6,1)=  Limiting  ply  shear  stress 

ALL0W(1,2)=  Liinitinq  ply  compressive  strain,  parallel  to  fibers 
ALL0W(2,2)=  Limiting  ply  compressive  strain,  transverse  to  fibers 
ALL0W(3,2)=  Limiting  ply  shear  strain 

ALL0W(4,2)=  Limiting  ply  compressive  stress,  parallel  to  fibers 
ALL0W(5,2)=  Limiting  ply  compressive  stress,  transverse  to  fibers 
ALL0W(6,2)=  Limiting  ply  shear  stress 
ACC  = Iteration  accuracy  factor 

NEM  = Element  number  of  element  chosen  for  stress  and  strain 

outputs 
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CALL  LLFAIL 
oICP 
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*rkllL  (farbObO)  ((',Y(Ir.J)r  I = lr2),  J = IrNYF) 
lAO  r.HlTt-  (orbObU)  O-blMU*  I " l.'JEC) 
i.itlTt,  (OrbOoO)  C-lr  b2r  Ur  V12 
r.FIlL  (OrbCrbu)  ( r.ikuLr.  ( I ) ' I -■  IrNANO) 

..KllL  (OroOOlj)  (iMK(I)r  I = IrNAf.G) 

<.KITc  (UrCOlO)  ( (ALLO.v  (IrU)  r J = 1,2),  I = 1,0) 

Ai<  1 TL  ( o r b(j  Vo  ) ACC 
r'Klli.  (G,bObO)  IJLM 
r'H 1 1 u (6,6020) 

"C.TUKM 

1010  FO'If-'/yT  (20  ( Ib,  lA)  ) 

1020  FOKMAT  (6 (bib, IX) ) 


I 


72 


•r  /I,  \ ? .*  ! \ ^ 


lOJU 

1040 

10  JO 
lUoO 
llj  70 
1 o6o 
lOOO 
5000 


f-CKMA) 
f-  ^ M'1  A I 
FOKMaI 

F 0 K I'-'i  A r 
t-  OkM  A 1 
i-^AI-AT 
F 1 
FoF’Ma  f 


bOlO  FOKMaI 
5u2o  FCaMaT 

1 


120F4. 1) 

( li'XALl'i.O) 

( 1CF6.U) 

(lOFO.l) 

(lOFl). 

( 1 OF U.b) 

(oFU.b) 

( /»  2X.  l^HFO'A  LFfcre.O 
1 OiiCHAUGLU  F(ji'.f-'  r I2i 
(iia,  2a, 

( //,  2X, 


T I-J,  24h,  The  OAfJOAlOTH  H/.b  bcLlJ  > 

4h  to  , 12) 

0 { 1 4 AI lE  r>T!!b5S")) 

>»‘*hFCK  l<^FE^'t^CE,  Tr(E  ItJPUT  DECK  I3  FEFKODC'CEO 


ITS  LFflRtl  r , //,  4X>  16A5) 


50oU 

F OKF'aT 

(lA, 

2016) 

5U4U 

F CKI-‘i  A T 

(6(5a,  314)) 

5 0 5 0 

(■  OttMAT 

( 6 ( 1 X 

' 2F7,4)) 

5Uo0 

)•  0 A f'i  A T 

( 15X, 

4E15 .6) 

5070 

FOAMaT 

(IX, 

1 C F 1 2 , o ) 

5 u h 0 

foknat 

(IX, 

10F12.1) 

5090 

F uKi  1a T 

( 5X  , 

lUFa.5) 

bOOO 

F Oki’T  A I 

( OX  , 

lOFfl.b) 

oulO 

FufU-lA  T 

( 6(5 

A , 2i_  1 4 , o , / 

DOtO 

F OKi'IAT 
ENO 

( //2X9  Uhi-i-PlANl 

» ) 

'oTFEbS-)  ) 


tine  S_t^i_f g n 


boiiKoUTiue:  STiFGtj 

coMMOH  / cor,sn  / 

CON. .’■.OH  / ..ATUAl  / 

CoN'Mulj  / ANCjOT  / 

CCMMOti  / CONiMl  / A^o-:F(o) 

v,;l  = VI r.  ♦ l2  / tl 


f',F  , 
Ell' 
(-.1.  S 


i^EX, 
L < , 
C.  ' C' 


(jSYf  HAf'Oi  HFAII.'  HXl,  f)X2,  fJTl,  NY2 

(.15'  E2"'..  F25,  i.N.J,  EE  (5, 5) 

V.12,  ri'Mb),  AI,:LHb)<  ALLOa  (6,2)i  HT 

.(5),  iiraso'b) 


0(1,1) 
0(1,2) 
o ( 1 , J ) 
0(2,1) 
0(2,2) 
0(2,5) 

0(5,1) 
0(5,2) 
j ( 5 , 5 ) 

nr 

Ui 


El  / ( L .0  - Vl2 
V 2 I ♦ C ( 1 , I ) 

(E2  *V1,:.*V21. ) / 

o(i,2) 

L,-.  / (1.0  - /12 
(U  ♦Vi,:-tV2i)  / 
0,0 
0,0 

2,0  ♦ G 


T 0. 

♦ V21) 


(1.0  “ 

♦ V21) 
(1.0  - 


V12AV21) 


V12+V21) 


+ 2.  * G(l,2)  + 2.  * 0(5,3)) 


0, 

.125  * (5.  ♦ (0(1,1)  + (J(2,2)) 

02  = .500  * (0(1,1)  - 0(2,2) ) 

U5  = .125  A (u(l,l)  + 0(^.,2)  - 2.  * 

04  .125  ♦ (C(l,i)  + 0 ( ,. , 2 ) + o.  a 

U5  = .125  * (5(1,1)  ♦ C(,.  ,2)  - 2,  > 

uO  10  1 = l,riA()o 

AN02i((I)  b AHOLud)  ♦ ((2,0  ♦ 5.141592653)  / 160.0) 


0(1,2) 
0(1,2) 
(V  ( 1 , 2 ) 


“ e.  t 

- 2. 

4 2. 


0(5,5)) 
0(5,3) ) 
0(5,3) ) 


10  Ml  : NT  f ImK(I) 

t mE  IFAHbF OKi-iED  uAMIhA 
Ell  = 0. 
bl2  = 0. 
cl5  = 0. 

C.22  = 0, 


stiffness  N'aTRIX  (E)  15  COMPUTED 


12 


pt,  T * 51  n>?  iT 


= 0. 

L.JJ  = 0. 

DC  20  I = lir.At,^ 

uli  = Lil  + (U1  +U2  Ar:02K  ( i ) ) fU,^*C>^S(2  .-*Ar.02iU  I ) ) ) ■•  IMI',  ( I )/hT 

cl'd  = E12  + ( U4-UJ)  > l.ui>  I .»  A,  |02K  ( I )))♦  THK  ( 1 ) /MT 

h.lO  r Li-5  + (0  ,bO*U2*'SI,.(  AL02K  ( 1 ) ) tbS^SI'n?,  ^AI.G^K(  1 ) ) ) ♦ Ttlr- ( I ) /rT 

l22  = cl^2  + ( Jl-.^L  >COS  (mI.G2-<  ( I ) ) ♦U-5‘i:OH2.*A  a;I.in  I ) ) ) * lit*'- ( 1 )/llT 

l2-5  z tii-J  » (0.‘ju<l£  t'.ai.l  A'i-.-c:iM  1 ) )-03*LIL(2. 'AI.02K(l  ) ) ) ( 1) /liT 

u-J-S  = LOO  f (Ub-U.i*i.Ob  , iAI.02iM  I ) ) ) ♦Tllr<-(  n/nr 
20  COMifJUt 

rtHlTt  (o.bOOO) 

nuiTt  (u(bOIO)  tllf  Ll2t  Ll3r  t22r  F2i»  L33 
KLlOKf< 

bOOU  LOkN.aT  UHli  IbX  tbUlMt  lOMPOSITL  At'ATRIX  IS) 
bUiO  LOKMAI  (/IbX-iLlb.b  //  00A2Llb.b  //  AbXElb.b) 
lGl) 


Subroutine  Constr 


aUbKuUTiuL  CONSTR 

C0NSTi<  FIRST  CA.i.S  ELCC.JS  V'lICU  C OrjS  1 RUC T5  ThC  r.OIvIOUAL  rLL"LMT 
jlI)ri-L.SS  PA  IR  iCi-.b  Ai  ij  Thl'J  U'.iLS  TmL  LLLLLM  F'AH  ICt-S  TO  CONSIHUcT 
I Mr.  ol  IFrM.bS  ' A I ('■  i KOi<  THS  'vHCLr.  STKUCTURfJ, 

Cor-’.'or.  / ARKAYb  / A(00.2u0),  FO(20r))f  ri')(?0U)»  L IS(200)  ,T[)lS(20u) 
COMPuh  / «.LC''-'LT  / IbJ>  Nl.L,  i.KDi  NZC.  M:'(4SC).  l-..'r”(  >0)»  XY:'(20C) 
lO  lu  J = lilIKO 
F 0 ( -J ) = 0,0 
UO  Ju  1 ItloO 

A i i » J ) = 0.0 

lu  CONfil.OL 

00  20  I = liNLL 
IS  2 I Y -J 
CmLL  LLCCrjS  (1.5) 

20  COilTir.UE 
)<L  T'JKN 
LNU 


Subroutine  Elcons 


SOoROUTliu-  Ll.COfIS  (IJ) 

CO“-;Mol.  / ARRAYS  / A (.50.200  ).  FO(200).  ri;;(2C0).  D I S ( 2 0 0 ) . TU 1 5 ( 20  0 ) 

CO;*.Mo;i  / GLvy.vCT  / ibU.  NLl.  . f.Rl).  IIZC.  nM(4'.5n.  l,/'p(40).  XYM(200) 

COiA'-'.ol < / ' ' A I u.T  r / L 1 .1  . t,  i 2 . L 1 -T I [■. 22  . L2 0 . f .0 0 . 0 ( S . .5  ) . lc.  ( -5 . .5  ) 

ui  -'-USlOr.  (.02  (-.5) 

I.Ou(l)  = l.f.;(U-2) 
z i.r.'(  I )-i ) 

i:ou(S)  = i.r-'(  u ) 

CO  20  I = 1.2 
M = I + 1 
(.0  20  j = M.3 

IR  (iJ0L(  I ) - t;00(  J)  ) 20.  10.  10 
10  (IT  = fiwl,  ( ( ) 

NUU(I)  = MCU(J) 
i^Oj(o)  = riT 
20  CONTINUE 
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f.OU(i)  * 2 

♦ 2 

NOUIJ)  ♦ 2 
NIY  - I 
rj2f  - 1 
NJr  - i 
u2r  - hiY 
NJl  - IJIY 
IJJY  - ,j2Y 
Xfr-.(NlY-l) 
XY-<(ra  f-i) 
XTM(N2f-l) 
XYMdJlY  ) 
XYM(f>ilY  ) 
XY.'(N2t  ) 
Aijb  1X12  * 


t • r i T.* 


r rnoN 

t Uit  ' 


) - XYM(fl2Y- 
) - XdvdJoY- 
) - XY-'dUY- 
) - xrt-‘dj2Y 
) - XYrddJY 
) - XYM(d.)Y 
YiJ  - Xlj  ♦ 


11 

1) 

1) 

) 

) 

) 

Y12)  ¥ 2.0 


r ( AK4  ) 
xi^;  / SAK 
X.IJ  / sak 
x2.)  / bAK 
Y12  / bAK 
YIA  / SAM 
Y23  / bAfi 


YIA  / 
Y23  / 
Y2A  t 

Ai,}  » 
Xt  A > 
YU  » 
XI A ♦ 
XI A A 
Y12 
X12  ► 

A 1 .i  ♦ 

YIA  » 
Xii  * 
XIA  ♦ 
YIA  ♦ 
Y12  * 
X12  • 
X12  ♦ 
Y12  ♦ 
Y12  ♦ 
X12  ♦ 
X 1 2 * 
Y12  * 

1 t M X ) 

2 t i . 1 A ) 

■ 1 n. 1 A ) 

2 f 1. 1 X ) 

• 1 <dlx ) 
•2*;. IX) 

1 il.lY  ) 
...ir ) 
1 n.i  r ) 
*,.1Y) 
►l»r,lY) 
1 ||,2X  ) 
2»d2X) 

^ 1 f ^ j 2 A ) 
► 2 ditX  ) 
1 d.tY  ) 

• I.2Y) 


SAM 
bAfi 
Y2A 
;<2A 
Y2A 
VIA 
XIA 
YIA 
Y 12 
X12 
Y12 
Y2A 
a2A 
Y2A 
A2A 
Y2A 
x2A 
Y2A 
X2A 
YIA 
XIA 
YIA 
XIA 
= A( 


A ( 

A d.21 1-1 
A ( ' j2  1 ¥c. 
A d < A U 1 
A (mAI  t .. 

A I 1 
A (d21 
A in.'  I U 
A d.Al 
A dUl  <-1 
A ( 1 

A I 2 
A dU2  t 1 
A d>iA2  t2 
A ( 1 

A (NA2 


+ E1UYY2A 
+ ElAtYY2A 
+ L1o*XAAY 
+ E.lcvXAAY 
4 C.114YY22 
4 tlA'Y'd;!; 

4 f.;r.4XX2A 

4 t, i.:‘yaax 

4 r<-A>YAAX 
4 E2A»;<X22 
4 F.22«XX22 
4 (1114YY1A 

t LIA.YYIA 
4 UA*X11Y- 
4 L12»X11Y- 
4 F.2t»XXlA 
4 r.i2*mx- 


- 2,'J1E1A«XY2A  4 lAA4aX2A 

- (E12»EAA) 4XY2A  4E2A4Aa2A 
L 1 1 * Y 1 0 A •-  L A A 4 A V A A 4 1 1 o + ir  A A X 
i;  134  Y t A i-L<  AvVXAA  I i-.Av  ‘ YAAX 
4 i-ji  '>,.'2  (■  -2  1 A ^ Y .'.id.  ( 2 A J ♦ i X22 
H 2 4 .U  Y - L . i A 4 V .1  :f  > 4 r.  2.  J 4 X X i <- 

- 2.0*i;2A»XY2A  4 LAA*YY2b 
1 2 A ♦ X X . j - L 1 b 4 Y Y .i . j 4 L A .)  > X A A Y 
r. 2 * X).  A A - 2 A A 4 V r A A 4 c 2 A 4 X A A Y 
L124  Y22X-i,AA*X2?y  4-E1A4  YY22 
r.2A*  Y22X-L2A4X22Y  4L  AA4Y  Y22 

- 2,041  1A*XY1A  4 lAA4XX1A 

- dA124EAA) 4X1  134l2A4aY1A 
■t  1 1 * Y V 1 1 -t  3 A 4X  X 1 1 44.  U 4 Y 1 1 X 
f.lA*YY  1 1-l234  vxll  41-..104  YllX 

- 2.0*L23*XY1A  4 lAA*YY1sj 
■C23*XXll-E 1A*Y Yll4tAA»Xl  1 Y 
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r>  — * '*  ? V 5 s>  ? r rODV 

.wfiii^iULs;  Lui  I 

/*  (r<5£>+l  ,:,•  Y ) r A (.jJ'.- +1  tu^v;  ^ Y 1 1 Y-c.,-J?  v ,.  ,■,  n 1 m < f ^ j5»  X U Y 

m(  = A(  i,t....X)  ♦ L11‘YY12  - IJifXYl^  t ,;.V5*^X12 

Al  2.I.JX)  = A(  2.t(jX)  + 'iLi5>YY12  - U:  1 2 ‘ c ) ♦ X r 1 2 -*  c.-' 5 • aX  1 2 

A(  ltl,.iY)  = Al  I.IIJY)  + L22»XX12  - 2 . 0 < ;,2  j i a Y 1 2 f L3J»YYJ12 

Ac  1 UKN 
LhO 


Subroutine  Excite 


bUi.KOUT  If,t.  EXCITL 

COMMur.  / Ai<HAYS  / A(30»20C),  FO(20(1),  rin(200).  DIS(200)  iTD1S{2.0(j) 
CJf-Y-iOts  / oLO;-LT  / UUi>  lJ2Cf  tw>’('*0).  XV-'(200) 

COi-iMON  / aLOAuS  / r.Xitr  mXF,  (,br’X(10)»  NrP;<(Uj)t  l.X(2flO)t  f;A(2) 
C0>:I'10I-  / YI.OADS  / ti'CL),  l(YF-,  ij:;PYtlO)-  NFPY(10)»  i,Y(2,10)i  EY(2) 

Idl  = Ibu  - 1 
hCjC  = nxu  + UYO 
iF  (NU'C  ,LU.  0)  uO  TO  2HU 
u 0 i 3 0 I = 1 I f . U'  C 
IF  ( i - I.Xtl)  iU»  10,  20 
10  Jr.  r.LHX(I) 

u = t-xm 

Jo  iO  JO 

20  J = NOPY  ( I - NXO) 

0 = LT(l) 

JO  IF  (IbJ  - J)  40,  'tO,  60 
40  DO  50  ^ 1 , Ibi 

FO(j  + K-IUD)  = (■vjlJtK-U-t;)  - A(l“K  + IUO,JfK-I!iD)  * U 
bo  cciiririuc. 
oo  ro  JO 

60  IF  (J  ,tO,  1)  Go  TO  60 
JM  = J - 1 
jO  70  K = 1,JM 


FO(K)  A FOU)  - Al 

vJ-K  + l#K)  * D 

/o 

cowlifiut 

60 

IF  (KUJ  - 

rjAU  f o) 

90, 

90,  110 

90 

JO  100  A 

= 2 , I L)  iv 

FOl JtA-1) 

r F G t J + K 

-I) 

- A (K, J) 

lOO 

cur,  r IM'::. 

uO  TO  .100 

110 

IF  (J  ,t_0. 

IjRti)  GO 

To 

IJO 

I L z N K J “ 

J *■  1 

GO  120  K = 2,  IL 

FO(vtX-l)  = F0(JH<.-1)  - A(K,J)  X D 


ic-.O 

COhr  Ji,G2 

IJO 

COPT IKUt 

GO  2JC  1 = 1,NDG 
IF  (I  - i.XU)  140, 

i'*G, 

ISO 

loO 

,j  = flCPA  ( I ) 
D = EX(i) 

GO  ro  160 

IbO 

J = hU'Yli  - NXij) 
J = tYll) 

1 oO 

IF  (160  - J)  17u, 

1 70, 

190 

170 

DO  IbO  X = 1 , It.l 

A(  1-X  H,;0,«.,  + K-IrjU)  = 0 

.0 

160 

CCIlTINUL 
GO  TO  210 

190 

IF  ( J • t-0  , 1 ) 00 

TO  210 

r 


[ 

I 


/f> 

rv/Y 

oM  = J “ 1 
00  200  K = liJM 

A I vj-K  ^ I I K ) = 0,0 

2UG  CONflr^Ut. 

21U  bO  220  K z 2(idU 

A(^l  J)  0.0 

t2o  courujuL 

K0(J)  = All.J)  * 0 
,.00  corjTii.uc 
2>+0  in;I-C  = NXf-  + NYI- 

IF  (bFC  .to,  0)  bu  TO  2o0 
00  270  i = I.bFO 
4F  li  - t.Xr  ) 2!>U.  2b0*  2bQ 
2b0  ..i  = f\FFA  ( I ) 

K - l'iX(iil)  ♦ 2 

L = NX(2t  I)  ♦ 2 

L)  = FXU) 

b - (XfM(L)  - XTM(K) ) / 2,0 
bO  ro  270 

2b0  j = I.FPT  ( I " HXF  ) 

K - (1X1,1  ••  IIXF)  ♦ 2 - 1 

L = NY  (2, 1 - l<IXF  ) ♦ 2 - 1 

0 = LY(2) 

S = (XYM(L)  - XrN:(K)  ) / 2,0 
t:70  FO(j)  :i  l-'O(j)  + b + D 
260  CONTlIiUii 

CO  oAO  1 = i,N2C 
J = N2PU) 

FO(J)  = 0,0 
00  ^'JO  K = 2,IHU 
A(K,J)  = 0,0 
,;90  COiUINUt; 

IF  (IbO  - J)  000,  300,  320 
300  UO  310  K = l,Ibl 

<,  ( 1 “X  > 1 !;IJ  , ,J4-K  — I 60  ) " 0,0 
31C  CONTKjUl 
oO  TO  3'*0 

320  IF  (J  ,c.O,  1)  Gb  TO  340 
JM  = J - 1 
bO  330  K - 1 , JFi 
330  A(1-K+J,K)  = 0,0 

040  COfNiIl.Uc. 

bO  3b0  1 = l,l,Hb 
F1N( n = FO( I ) 
obb  continue 
Kfc.TU,(N 
2Nl) 


^rubrou t i ne  Cause  1 


bUbHOUTI.'.u  OAJbCL 

COMMUf*  / AKKAT'.>  / A(JO»aoO)»  P0(<;00)-  FlfJ(200)f  u I S ( 20  0 ) . TD I b ( 200  ) 
COMMuf.  / oLO.vLf  / lU.^.  i<LL»  ta(0»  N2C  > NV.(4bO)r  !,tM(4U)»  <YM(?.00) 
lut.  = Ibu  “ 2 
la  = I.KU  - 102 
r,2  = ra  + i 

bb  au  N = 2 1 M l 

J2  = 102  + N 
UO  10  I ::  N.J2 
K = AU-ta2,!j-l)  / A(1>M-1) 

A(I-N  + i.M)  = Ad-Mi-lffa  - /<  ♦ A(2»M~1) 

(-0(1)  = F0(I)  - f<  ♦ FO(i;-.l. ) 

10  coraiMUt 

M = M -*•  1 
UO  20  J = M,J2 
UO  20  I = J.J2 

A(I-J+l,J)  = A(I-J  + 1»m)  - A(I-M  + 2,M-1)i.A(J-M  + 2»M-1)  / AU/M-l) 
20  oOMTlMUc. 

UO  HO  M = M2,rJKo 
UO  30  1 r N,J2 
r<  = A(  1-Mf2.;j-L)  / A(l.M-l) 

A(I-MHl.M)  = AU-Md-M)  - h * A(2<M~1) 

KO(  1 ) = FOm  - K * F0(N'*1) 

30  comiiijuf: 

M ::  M •>•  1 
00  HO  J = M,J2 
UO  40  I = JtJ2 
IK  (J-i\RD)  45i  4!)»  35 
35  A( I-jHl, J)  = 0,0 

45  A(l-o+l|J)  = A(.t-J  + l*J)  - A(I-M  + 2pr:-l)  ►A(J~)J  + 2»!J-1)  / A(1,M-1) 
40  UOtait-UL 

UO  oO  J = lillRQ 

H = 0,0 

I = f.RQ  - J 
00  50  K = 2,100 
IK  (I+(\-InRO)  55,  55,  80 
80  TOIS(Ihk)  = 0,0 
Db  R = i\  f AIK,I  + 1)  + TOIS(IhK) 

50  COi'iriMUL 

TUIS(1  + 1)  = (FO(iHi)  - (n)  / A(1,I  + 1) 
oo  CONTiMUt. 
kLTURM 
LNO 


1 


Su_brou  t i ne__E  l_f  a ij 


bJbKOUTLME  Ll.KAIL 

LONlMor^  / AKI<AY5  / Al30.2U0)»  F0{?00)»  FirM?00)»  IJI5(2C0)  ,TOIb(200) 

CO^'.MO^  / CCUbfl  / ritb,  Iibx,  f.bYi  :!A'.5.  r;FAlL;  I.Xlr  U>.2f  I.YI,  fiY2 

cuMMor;  / bto.-Lr  / ibL'f  f.'bUf  f.RD»  .‘j/c  i r,/:P(‘*o)p  xy:-'(200) 

CbPil'iOfp  / ‘•MT-.'Af  / tu»  1.12»  Flip  f.22p  L2b,  F.ob,  4(bpb).  I.FC-ipi) 

LOK-;uf,  / Alir^OT  / f.lp  r.2p  Gp  V12,-  ;HK(b)p  Ar;;.Lb(b)i  AbLO,,  ( 6 p 2 ) p MT 
CO’-,N'um  / n / AHo2f<(b),  rHFUCbjp  ITT(lbO,b) 

CGHi-lorJ  / CC-t'T2  / AC^f  NlM 
ulf.LliSlO,,  5ri!t.(b),  brl^S(G),  F00(200) 

UO  4 1 = IpilANG 

THlIA(I)  = AbOLL(I)  ♦ (<5,14159  / lOU.O) 

4 COlpTlI.UE 
IP-pC  " 1 
b OONTlriUn 

(jO  10  I = IpNRQ 
FOU  ) - 0,0 
10  COtjIll'JUb 

uO  00  I = IpULU 
10  = I ♦ 0 

CALL  SIRAIli  (10,  LPX,  Ef'Y,  OXY) 

CALL  LAfEH  (10,  tPx,  LtY,  OXY,  LAY) 

IF  (LAt  ,Nl..  0)  call  PFORCE  (IO,  lPX,  EPY,  GXY) 

IF  (1  .f-L,  ipEM)  bO  ro  OU 
LL(1,1)  = Lli  - cL(l,l) 

£.1.(1, 2)  - Ll<i  ~ Lc.(1,2) 
lL(1,0)  = F.lA  - EL  (1,0) 

LL(2p2)  n E22  - c.E(2,2) 

t.t(2,0)  = 1,20  - tt(2,0) 

£-L.(OpO)  — c.0a  - Ec(0,0) 

t.t;(t:,l)  = tt.  (1,2) 

Ll(0,1)  = Lt(l,0) 

t-L(0,2)  = Ei-(2,0) 

UO  2b  K = 1 , 0 
bfp<b(K)  = 0,0 

20  continuf: 

blKN(l)  = LPX 
bIRN(2)  = LPY 
bTRN(O)  = GaY 
LO  27  K = 1,0 
CO  27  L ::  1,0 

bTKS(^)  = STHb(K)  + Lt(K,L)  * bTRN(U 
27  COItlKpbt 

00  20  K = 1 , 0 
Go  2o  L = 1,0 
t-L  ( K , L ) - 0,0 
26  CwIpTINUL 
00  CONTINUE 

uO  Ob  I : 1,U2C 
J = li2P(  I ) 

FO(J)  = 0,0 
Ob  CONTK.Ul 
Fo'''Aa  = (j.O 
GO  40  I : l.NRD 
F = F0(  I) 

F0(i)  = FO(I)  - FOO(I) 

FOO(I)  = F 


FOMaX  - AMAXl  ( Au!>(FU  ( 1 ) ) > FOMAX) 

4 0 COfdlNL'L 

.SilTL  (o,4'-j)  IMC.  IITr  FuMAX 

4b  FoKMaT  (U.Xi  UH1..CKl.  ■'MiT  =f  Ut  1IU»  2Xf  1 211 1 f [[.RAT  I ON  =r  ilii 
ilfirf  2X»  NIFOMAX  r:  r ti  i 1 • A ) 

:f  ifonax  ,lt.  ao^!  uo  to  so 
in  = II  f + 1 

IF  (I  IT  ,0f.  20)  bO  ro  70 
CALL  FORr.UU 
bO  TO  b 

bu  calL  output  1 ir;c»  nr) 

.UaTL  (onS)  NLM 

■lIvlTE  ( L)  I r 0 ) (bllirj(l)»  I - 1(3) 

mRITi-  (orhb)  NtM 

..rUlL  (O(<’0)  (bTK'o(I)/  I = 1(3) 

.'.KITL  (odUO) 

J-i  ~ I'JhU  / 2 
UO  Ob  I z 1(M 
J = 2 * I - 1 
X = 2 ♦ I 

rtKITL  (0(110)  I(  TU1..(J)(  TuI5(K)(  F0(J)(  F0(I<) 
bb  CONTiNUL 
v.RlTL  (0(99) 

INC  = INC  ♦ 1 
00  OO  I z 1(NHD 
F0(1)  = FIN(I) 

00  COi-JT  INUl. 

CALL  FOKr.DU 
IIT  = 0 
oO  TO  b 

70  CALL  OUrpUT  (INC(  IIT) 
mKITc.  (6(7b)  NrlM 
(,SITL  (U(O0)  (SThU(I)(  I = 1(3) 
i'iKITL  (b(Ob)  NLM 
AKilc  (0(90)  (STHb(l) ( I = 1(3) 
oRITL  (odOO) 

Fi  z nKD  / 2 
CO  72  I z 1(M 
0 = 2 ♦ 1 - 1 
K = 2 ♦ I 

whirc  (odlO)  1(  TDIS(J)(  Tu1S(K)f  F0(J)(  FO(K) 

7c  cOfnif.UL 

7b  FoRNAl  (lhO(  lbX(  cF.mM'L  STRAIN  IN  ELEMENT  ( 13(  2X(  SHISd  //) 

oO  rOKi-i,vl  (1/X(  10)(X-bTi<Al.<  =(  L12.b(  / 17Xi  ICnY-bTfiAir-  Z(  El2.5( 

1/  17X(  l4HSI!LAI<  Sn<MlrJ  =(  E12.5) 

Ob  format  (1mO(  ibX(  22Hri'i;  STRESS  IN  ELEMENT  ( I3(  2X(  3U1S:(  //) 

90  FwKMhI  (d’X(  lOHX-STivuSS  Z(  tl2.b(  / 17X(  lOMY-STRESS  Z(  F.n,5( 

1/  17X(  I'ihSmEAK  SIRESS  =(  EI2,b) 
l6o  Format  (i^^0(  lbX(  4h-,CDi.(  4X(  6)(X-DISP(  4X(  oHY-DISP(  5X( 
1/HA“FCRCE(  0X(  7hl“F0KCE(  /) 

99  F0(<>1AT  (IHl) 

110  format  (16X(  13,  2(4X,  F0,4).  2(4X,  Ell, 4)) 
kL  rURN 
EM) 


i 

i 

I 


80 


V i 


at  copv 


Subroutine  Strain 


bUbHUUTlIiK  STRAlIvi  llj.  Ef^X,  LPX,  GXY) 

uOKMOr.  / AKtOiYb  / A ( 30  i .;O0  ) t fC{200).  rr-'.OOO).  0 I b ( 200  ) , TO  I j ( 20  0 ) 

COMMOli  / utO.'-'EI  / lbu'»  iJLL>  rif<U»  .‘IZCi  fjV('t50)»  f;ZPl‘40)»  XY'-‘(200) 

laV  = l.i-U  13-2)  f 

liZY  = ♦ 2 

I.3Y  = 0,'-Ul3  ) ♦ 2 

X12  = Xff-.(NlY-l)  - XYI',(M2Y-1  ) 

Alo  = XVi-iraT-i)  - ,\fO(i)3Y-l) 

a23  = XY;’(M2t-i)  - xr:i(,j3Y-l) 

Y12  = xYMuar  ) - xi;aiM2Y  ) 

Y13  = XYnUar  ) - xr'(f(3Y  ) 

Y23  : XYO  (W2Y  ) - X'iV,(N3Y  ) 

nR2  = ALlblX12tY13  - XISaYIZ) 


Ui  =lElb(iaY-l) 
u2  =TLlblYaY  ) 

03  =lLab(u2Y-l) 

=I0IS(!i2Y  ) 

UU  =rOIb(l.3Y-l) 

L o = 1 u 1 J ( 0 3 Y ) 

EHX  = {Y23»D1  - Y 13*03  + Y12*0‘3)  / AU2 
LPY  :-(X23»02  - Xi3*L‘t  * X12*Ob)  / A!<2 

oxY  = (Y23  + 02  - Yi3»uA  + Y).2<0b  " X23  + D1  + X13»D3  - X12*U5)  / AK2 
KtrURf) 

LOO 


Suty^outXfL®  Layer 


bULlKoUTIOL  LAYER  (I3,  Lt'X,  LPY.  GXY*  LAY) 

LOV.'-luri  / COfJbTl  / OL3*  Ubx,  OSY,  *'A':G*  O'^AIL*  TJaI,  fa,',  ,aa.  MY2 
uur-MoO  / .uvrUAr  / LU*  L12,  L13.  C22*  r.23,  E33,  C(3.3).  tL(3,3) 

/ AfaboT  / Ll<  E2*  G*  V12,  THn(5)*  AnGLE(b).  ALLG.aOiZ)*  )|T 
tOV.MOi!  / CC'VfMl  / A).u2R(-j),  IrETA(b)*  111(100*0) 

*<LAL  c.blKO(3),  LI-.brRo(3),  LfbIRb(3),  T(3*3) 

tblRiNd)  - lPX 

t.brRi*(2)  = LPf 

i.bTR,'a3)  = OXY  / 2.0 

LAY  = 0 

0=  I3  / 3 


uO  1‘jl 

- l*IJA(,o 

ITT ( J * 

K) 

= 1 

1 d*  1) 

:: 

(COS(Trlc.TA(x)  ) ) ♦ ( COS  ( T)|F.  T A ( K ) ) ) 

1 d*.i) 

= 

(b:(.(lM(dA(K)  ) ) * (SII;(ThFTA(K)  ) ) 

T (1.3) 

2,0  * bira  rhET*.  (K)  ) ♦ oOSdMETA(K)) 

T(2,l) 

:: 

T d * 2 ) 

1 (2,c:) 

- 

T (1 , 1 ) 

r (2,3) 

z 

- Id*3) 

1(3,1) 

z 

- bI:idMf.TA(K)  ) ♦ C0S(THLTA(K)  ) 

1(3,2) 

z 

- 1(3,1) 

] (3,3) 

z 

1(1,1)  - T(l*2) 

uO  10 

1 

: 1,3 

LMblR.J(l)  = 0,0 
Oo  10  L = 1*3 

LObThial)  = LI'bTR.Jd)  + rd*L)  » ESTKN(L) 
lU  OOOdnUL 


•} 

i 


8) 


(jO  t’o  I : l.a 

Lf-.jIK'ot  [ ) 0.0 

uu  <;u  L - 1..'! 

Lr-'.'oIHS(  I ) = LN'blRS(I)  + O(I.L)  ♦ Lf-'STRrJ(U 
d.ij  ocruifJUL 

L>ibIfvS(J)  = ♦ LMbTRM(O) 

DO  I'H  1 = 

IH  (Nf-/«1L  .LO.  ■/.)  00  10  bO 
If-  (L,'10TKfi(  1 ) .1,1,  0.0)  GO  TO  40 

IF  (Ll-iOTKN(I)  - ALLOV.(i,D)  140  , 140  , VO 

40  IF  (LI-ijIKI,(  I ) t ALLoW(1,2))  70,  140,  140 

bO  COiwT  IilUE 
L = 1 4 J 

IF  (LMbTKS(I)  .lT,  0.0)  GO  TO  bO 
IF  (LS'.bTKS(I)  -•  ,vL.uOW(l,  1)  ) 140,  140,  /O 

oO  IF  ^L^bIt<b(I)  + ALL0v>lL,2))  70,  140.  140 

70  CO;j(INuI. 

GO  10  (GO,  90,  100),  I 
00  IT  I ( J,K)  =2 
GO  TO  140 

90  IF  ( III (J,K)  .Lo,  1)  Go  TO  9b 
I r I ( J,K)  = b 
GO  TO  140 
9b  ITI ( J,K)  =3 

oO  To  140 
100  M r ITT(J,K) 

oO  TO  ilOb,  lOo,  107,  140,  100),  M 
10b  ITT ( J,K)  = 4 

GO  TO  140 

106  JTMJ.M  = 6 

GO  TO  140 

107  1TI(J,a)  = 7 

GO  TO  140 

lOU  lTT(J,r\)  = 8 

140  COiJilNUu 

141  CCIITINUL 

IF  (1TT(J,K)  ,Lvj.  1)G0  TO  IbO 
f)  = ITT  (u,K) 

LAY  = LAY  + 1 
IF  (LAY  ,GT.  1)  GO  TO  14b 
OO  142  I ::  l,J 
00  14,;  t.  = l,J 
ti_(I,L)  = 0,0 

142  COlilli'luL 
l4b  COlTTINUt. 

CaLl  DlGTIF  (N,K) 

IbO  COJ.Tir.Ua 
Ibl  ^.uliruiUE 

lE(^,1)  = Lf-(1,2) 


Subroutine  Dlstif 


i. 

i { 


■jOuKOOrir.L  ULSTIK  (U.K) 


B 

CU^:MO,.  / V 

AI 

oA(  / f: 

11,  li2,  i lb, 

, E2b, 

Ebb, 

r 

(b,b) , EE(b,b) 

B 

OU''MoiN  / A 

I ti  'jv'T  / “ 

1,  L2,  , V12,  THaC;),  A 

:.oLi.  ( 

) , LL 0 ,1  ( 6 , 4 ) , ) (T 

B 

cor.N'o:,  / cor- 

on  / A 

(nr-Mo),  Hit  lA(b ) 

, irrabo,b) 

B 

0 0 To  (IOC 

f 

10,  0 1 

,-4  , b.'  1 bO  • 20  ( 

bU),  0 

B 

10 

>j  1 = , 1 b 

lb.  * ( 

0(1,1)  ♦ C.  ( 1 , b ) ) 

» 2.  ♦ 

0 ( 1 , .2 

) 

+ 2 . 4 0 ( b , b ) ) 

■ 

0'2  ~ ,bOo 

♦ 

('.  U,  , 1 ) 

- 0 ( i , b ) ) 

B 

UO  = .120 

(4(1,1) 

+ 0(1, b)  — 2,  ♦ 

0(1,2.) 

- 2. 

0 ( b , b ) ) 

H 

04  = ,l..b 

(0(1,1) 

♦ C ( 1 , b ) ♦ 4 , * 

0(1,2) 

— 2 , 

0 ( .( , b ) ) 

u 

0 -3  - , 1 2 b 

(4(1,1) 

, 0 ( 1 , b ) - 2 . ♦ 

0(1,2) 

+ 2 . 

4 

0lb,3)  ) 

K 

CjU  10  40 

B 

20 

ux  = ,i,;b 

♦ 

( 5 , , ((.(,•:,  i ) X y ( 2 , 2 ) ) 

4 2.  ♦ 

0(1,,. 

) 

4 2.  * „(s.,b)) 

B 

02  = ,bOJ 

4 

(0(2,0) 

" 0 ( ) 

B 

Uj  = ,12o 

4 

(-.(,:,  b ) 

i Q ( t-'  f . ♦ 

0(1,2) 

- 2, 

♦ 

0 ( b , b ) ) 

B 

U4  = .l2b 

♦ 

(0(2,0) 

+ t),  * 

0(1,2) 

-•  2 , 

* 

0(b,2)  ) 

■ 

o b = ,12b 

* 

( 0 ( 2 , b ) 

0 ( <- » 2 ) " 2 , ♦ 

0(1,2) 

4 2. 

♦ 

0 { b • b ) ) 

i 

0,0  To  40 

I 

00 

01  = .12b 

t 

( -J , , ( 4 1 1 , n -X  0(2,2)) 

+ 2 , * 

0(1,2) 

4 2 . * 4 ( b , b ) ) 

i 

02  = .bOO 

♦ 

(0(1,1) 

- 0(cr2) ) 

1 

ob  = .12b 

* 

1,11.1) 

+ C(c:i2)  “ 2.,  ♦ 

0(1,2) 

••  2. 

* 

0 ( b , 3 ) ) 

1 

U4  - , l2b 

♦ 

(.,(1,1) 

t-  0(c»2)  ♦ 0,  ♦ 

0(1,2) 

- 2. 

♦ 

0 ( b , b ) ) 

1 

0 D ::  .12b 

4 

(0(1,1) 

► 0(c'(2)  ■"  * 

0(1,2) 

4 2 , 

f 

0(b,3)  ) 

40  001,1  ir.ot. 

Lt  ( 1 ( 1 ) =t.r  ( 1 , 1 ) + ( 0 1 mJc  «CoS  ( Ai  K ) ) + L'3*C0S  { 2 , ‘Al { K ) ) ) * TmK  ( K ) 

1/Hf 

( I , 2;  M 04-,.  ( ? , • AIJ62R  (r  ) ) ) *TMK  (K  ) /mT 

C.L1  1,0,  + (0,oO,vJc:»Sl!jtANOcR(K}  ) ) ) * ItiK  (K) 

1/hT 

cL(2,2)=LE(,i.2)  + (ul-U2-tC0L-(,\'iG2K(K)  ) +U3*  C05  ( 2 , * Ar32K  ( K ) ) ) »TliK(K) 
1/hT 

01.  (2>^)  -Lt  (2i  ;j)  Mii.b  J •U2*SI!.  ( ANOcR  (K)  )-US*Si:.(2  . > ' U02R  ( k ) ) ) ♦ THK  ( K ) 
1/Hl 

Ll.(o»0)=LE(J.  J)  MU5-uO*COS(2,  >AI,02R(K)  ) )+THK(K)/Hr 

loo  Cull r ifiuc. 
htruKH 
UJU 


i 


f 

I 


A0-A0»3  7*9 
UNCLASSIFIED 
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air  force  INST  OF  TECH  WRiOHT-PATTERSON  AFB  OHIO  F/8  ll/4 

PROGRESSIVE  failure  OF  ADVANCED  COMPOSITE  LAMINATES  USING  THE  F~ETC(U) 
MAR  76  G E BROWN 


AFlT-CI-77-59 


NL 


END 


9-77 


Subroutine  Pforce 


' L.Ubf<OUTiriE  PK'Rcu  i-PX,  t:PY.  GXY) 

t COMMON  / Ar(PAY'‘>  / AU'IO.COC),  fO(-iOO).  FIN(2CU).  ulS  ( 200 ) , TOIS  (200  ) 

i COMMOtJ  / uPCi  £r  / Ib'L/*  l.-l.i  |JPD>  N?C»  tJ/P(O0)*  XYM(200) 

' COMMOfi  / 'lAIOAl  / LIU  tl3.  L22,  f_33,  0(3.3)*  LK(3,3) 

COMMOfJ  / /UII30T  / Ll.  L2.  G.  Vl2.  rHK(5).  AfiOLf  (b)  . ALLO,.  ( u . 2 ) i t(T 
1 l)  1i-)CN'->1  Cu  U ( o I 0 ) . bb  ( . 3 ) I LP(3) 

> CO  10  X - 1.6 

UO  10  J 1.3 
bti  ( I . J ) r.  0.0 
U (1.0)  = 0,0 
!(.  CONTlI.Ut; 

nlY  = NM(l3-2)  ->  2 
im2Y  = N.-K13-1)  ♦ 2 
I-3Y  = KM  (13  ) * ^ 

U(b.3)  = /.Yt-,(N2Y-1)  XYM(tllY-l) 

ij(3.3)  = XYMUaY-l)  “ XI,'-  (M.5Y-1 ) 

0(1. J)  = XYM(t,3Y-l)  - Xr,-(N2Y-U 

C(5.1)  = XYMdilY  ) - XY,.(fJ2Y  ) 

I 0(3.1)  = XYM(ti3Y  ) - XY,  (NlY  ) 

b(l.l)  = Xi;..(  .2Y  ) - XYi-.(N3Y  ) 

b(2.3)  - 11(1*1) 
c(2.2)  = b(1.3) 

U(A.3)  = b(3.1) 
ui(A.2)  = b(3.3) 

U ( b * 3 ) - b ( b . 1 ) 
c(6*2)  = (^(6.3) 
uO  20  1 - 1.6 
bU  20  A = 1.3 
UO  20  L = 1.3 

blKI.rv)  bo(i.K)  + l3(IiL)  * LE(L.K) 

20  CO.'jTii.Ub 

lP(1)  = LPX  ♦ HT  / 2.0 
t.P(2)  = cPY  ♦ HT  / 2.0 
LP(3)  = (.XY  * HT  / 2,0 
bO  30  1 = 1.3 

i KO(nlY-U  : KO(ult-l)  Y bb  ( 1 , 1)  ♦ EP  ( I ) 

I ^Od.lY  ) = F0(ulY  ) * i.'H(2.I)  ♦ EP(1) 

• Ku(U2Y-l)  r Ft;(i.2Y-l)  + iJj(3.I)  ♦ EP(I) 

P0(w2Y  ) = F0(.-2Y  ) + bU  (■*.!)  * EP  ( I ) 

F0(N3Y-1)  = F0(I,JY-1)  » L'BCj.I)  ♦ EPdJ 
KO(N3Y  ) = FO(iOY  ) + LH(b.I)  ♦ EP(I) 

30  CONTXNUl 
REIOKN 
ENU 


1 


84 


Subroutine  Foredu 


sudkoutii.e  fohluu 

COMKOli  / ARRAYL  / A(40t200),  FO(200).  rHU200).  DIS  ( 200  ) . TO  iS  ( 200  ) 
COMMon  / GtOi-l£T  / Ibu»  ncL,  NRU»  N4C»  tJM(450),  n2P(40)*  XYM(200) 
1U2  = lUU  - 2 
T4l  = i-IRU  - IU2 
N2  = la  + 1 
bO  10  fj  = 2.N1 
o2  = 1U2  + N 
uO  10  1 = N,J2 

FO(l)  = FO(l)  - (Aa-U*2,N-1)  / A(lfN-l>)  + FO(N-l) 

10  COUTIuUl 

00  20  i4  = N2tNRu 


(JO  20  1 = N,J2 

FO(I)  = F0(I)  - (A(l-r'H-2,N-l)  / A(1»M-D)  ♦ FO(N-l) 

20  CON  nr, UE 

UO  oO  J = 1»NRU 
K = 0,0 
I = NRu  - J 
UO  50  K = 2,1(30 
IF  (1+K.-M<U)  40,  40,  JO 


JO 

U1S(I+N)  = 0.0 

40 

R = K + A(K,  I + l)  ♦ 

UISII+K) 

50 

coMrinut. 

015(1+1)  = (FO(I+l) 

- R)  / A(1,I+1) 

60 

OONT  ii-iUE 

UO  70  I = l,HRO 

TUIb(I)  = TuIS(I)  + 

DI5(I) 

70 

CONTINUE 

return 

cNU 

bUb^OOTIiML  OUTPUT  (INC  I IIT) 

CU^•MuN  / CON'iri  / fcLir  NbX,  NSY*  NaNG*  NFAIL»  Na1»  NX2.  NY1»  MY2 

CuX.NiUN  / GbON'ET  / Il-ui  Nk.L.  IjPOi  UZC»  tJM(4bO)»  ii.lr'(40)»  XYN’,(200) 

COPrOI.  / aLCAUj  / N.<U»  NaF,  l.fPX(10)f  NKPX(lt)»  i1a(2»10)»  tlx{2) 

»,UNN(/N  / YLOAIjS  / riYU*  NYf  , (ji;PY(10)i  NFPY(10)»  i,Y(2.iO)»  F.Y(2) 

COF.F'Ot;  / AN  If, or  / Ei.  F2.  G»  V12,  THK(5)»  ANGLE  (b)»  ALL0.V(6,2)i  HT 
Cul-.NOr.  / lOMCU  / AN'j2H(G)i  1META(5)»  ITTUSOrb) 

Ic  (UT  ,GT.  20)  '.UatL  (u.lOO)  IIT,  INC 
IF  (IlT  ,LL.  20)  ..Hiru  (u#170)  INC,  IIT 
IF  (NV.X  .to.  1)  GO  ro  10 
IaLOmU  = tX(I)  ♦ FLOAT (INC) 

TILOaL  = LY(1)  » FLOAT(InC) 

•iMIt.  (6,110)  INC,  TaLOAL,  TYLOAD 
sjv>  To  20 

10  rxLOAu  = tX(2)  ♦ K.o,4r(tiiC) 

ncOAC»  = EY(2)  * FLOAT(INC) 

f.KHL  (6,120)  INC,  TalOALj,  TYLOAD 
2u  ..HITE  (6,130) 

•w\ itl.  (6,mo) 

if.Kili.  (6,lb0!  (AIJGLE(I),  I = 1,5) 

JU  JO  1 = l,NtL 

..KITE  (o,lfaO)  I,  (ITT(I,J),  J = l,i)MNG) 

30  CONTINUE 

100  F0,<MAT  (IHO,  15X,  29,itHE  LAMINATE  HAS  FAILED  AFTtK,  13,  / 15X, 

1 irm lEi'AriurjG  in  the,  i'»,  ix,  iohincpemeimT. ) 

110  FOHMaT  (1H  , 15X,  2JhJHt.  MAGNITUDE  OF  THE  APPLIED  / 15X, 

1 e4)tiJlbPL/.CE’'L.NT  IHKOUGH  THE,  I , / 15X,  13HINCREMLN1  WAS, 

1 F6.4,  o[U<.CHr.S,  / ITjX,  22HIH  THE  X DIRECTION  AND,  F6.4, 

1 / 15X,  26r<INCHES  IN  THE  Y QIHKCTION.) 

120  FGPMaT  (IH  , laX,  2bnrHE  'OAGNITUDE  OF  THE:  AF'PLIE.D  / 15X, 

1 ibHLOAo  Through  ihe,  i4,  ix,  ohincrement  / 15X,  3hwas, 

1 Fll.O,  IX,  14HLo/IN  uh  the  X / 15X,  OMFaCE  AND,  rll.O,  IX, 

1 12HLb/lN  UN  The  / IbX,  7nY  FACE.) 

130  Format  (IhO,  15X,  13,)FAILURE  CCOE:  / 17X,  14h1  = no  failure  / 17X, 
1 29H2  = failed  parallel  TO  FIBERS  / 17X,  34(^3  - FAILED  PERPI 

INDICULAH  To  FILEbb  / IVX,  19H4  = FAILED  IN  SHl.AK  / 17X,  47H6  FAI 
IllD  PARAuLEL  AuD  PERPINDiCULAR  TO  FIBERS  / 17X,  42,16  = FAILED  PARA 
li-LEL  TO  rlUERb  ANu  Ih  SHEAR  / I7x,  47tl7  = FailEu  FERPINUICDuAR  TO 
IFIDERS  AND  IN  SeilaR  / 17a,  29H6  = FAILED  IN  ACL  THREE  MODES  / 17X, 
1 2S/Hj,6,()  REFRuScNT  IOTAu  FAILURE  // ) 

14C  FORNaI  (22X,  4UHl«mI',A  1 LAMINA  2 LAMINA  3 LA'-’.INA  ')  LAMINA  5) 

laO  format  (15x,  7|-EuL”  > E(4HArjO=  , F4.0,  2X),  /) 

160  Format  (isx,  li,  /x.  s(  ii,  9X)) 

170  format  (lax,  IIHIHIS  IS  The,  I4,  IX,  lOHlNCREMEf.T.  / 15X, 

1 27hLDUiEloRIUM  w,.s  Oi)TAII,EO  IN,  13,  IX,  HHITERmTIONS,  /) 

RETURN 
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